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JiLl Abstract: We discuss extensions the CKKW-L and UMEPS tree- level matrix element 

^ ■ and parton shower merging approaches to next-to-leading order accuracy. 



The generalisation of CKKW-L is based on the NL^ scheme previously developed for e'''e~- 
annihilation, which is extended to also handle hadronic collisions by a careful treatment 
of parton densities. NL^ is further augmented to allow for more readily accessible NLO 
input. 

To allow for a more careful handling of merging scale dependencies we introduce an exten- 
sion of the UMEPS method. This approach, dubbed UNLOPS, does not inherit problematic 
features of CKKW-L, and thus allows for a theoretically more appealing definition of NLO 
merging. 

We have implemented both schemes in PYTHIA8, and present results for the merging of 
W- and Higgs-production events, where the zero- and one-jet contribution are corrected 
to next-to-leading order simultaneously, and higher jet multiplicities are described by tree- 
level matrix elements. The implementation of the procedure is completely general and 
can be used for higher jet multiplicities and other processes, subject to the availability 
of programs able to correctly generate the corresponding partonic states to leading and 
next-to-leading order accuracy. 
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1. Introduction 



Particle physics phenomenology has been awed by the accuracy of LHC analyses. The 
precision at which, for example, the Higgs candidate mass has been measured could only 
be achieved through a very detailed understanding of the structure of collision events in 
an environment that can safely be called messy. Remnants of single collisions alone give 
rise to large numbers of hadronic jets, leptons and photons - even before pile-up events are 
taken into account. In order to separate and determine the characteristics of rare signal 
events, highly accurate methods have evolved to describe background processes. 

Precise theoretical calculations for scatterings with multiple jets in particular are nec- 
essary for reliable background estimates. For generic processes this, until recently, meant 
that multi-jet tree-level matrix element and parton shower merging (MEPS) techniques 
were the method of choice, with CKKW-inspired prescriptions [1-5] being widely used. 
These methods impose a weight containing the parton shower resummation on tree-level- 
weighted n-parton phase space points. Phase space points with soft and/or coUinear par- 
tons in the matrix element (ME) event generation are removed by a jet-resolution cut, 
leaving only n-parton phase space points that contain exactly n resolved jets. The same 
cut is also used to restrict the parton shower (PS) to only produce unresolved partons 
as long as tree-level calculations for the resulting state are available. The combination 
of reweighting and phase space slicing (by the jet separation cut) allows to add tree-level 
samples with different jet multiplicity without introducing any phase space overlap. 

This method has an evident drawback, even if wc would be content with a tree-level 
prescription of multiple jets: Simply adding n-rcsolvcd-jet states cannot guarantee a stable 
inclusive (lowest-multiplicity) cross section. In particular, the inclusive cross will depend 
on the jet separation parameter, ^ms; so that choosing ^ms unwisely may result in significant 
cross section uncertainties. This problem is remedied by the UMEPS method [6], which 
infers the notion of parton shower unitarity to derive an add-subtract scheme to safeguard 
a fixed inclusive cross section. 

However, MEPS methods only improve the description of the shape of multi-jet ob- 
servables, and cannot describe overall normalisations or decrease theoretical uncertainties 
due to scale variations. This requires predictions of next-to-leading order (NLO) QCD. 
Through formidable efforts of the fixed-order community, we have recently witnessed an 
NLO revolution^, meaning that today, the automation of NLO calculations is practically 
a solved problem. Such calculations become directly comparable to LHC data by incor- 
porating the NLO results into General-Purpose event generators. NLO matrix element 
and parton shower matching methods like POWHEG [17-20] and MC@NLO [21-24] have 
- in parallel with the NLO revolution - become robust tools that allow a coalescence of 
resummation, low-scale effects and hadronisation with NLO calculations. 

The latest step in these developments are multi-jet NLO merging prescriptions [25— 
28]. These address the problem of simultaneously describing observables for any number 
of (additional) jets with NLO accuracy, and are thus direct successors of the tree-level 

^We cannot do justice to all results of these intricate calculations, so that we limit ourselves to the more 
conceptual papers [7-16], which made this progress possible. 
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schemes. The problem in NLO multi-jet merging is twofold. It is firstly mandatory - as 
in tree-level merging - to ensure that configurations with n hadronic jets are described by 
the n-jet ME. If we have a better calculation at hand, we do not want to predict rates for 
n hadronic jets by adding a parton shower emission to the {n — l)-jet NLO calculation. 
This problem has already been solved in tree-level merging methods. Secondly, each n-jet 
observable has to be described with NLO accuracy (if an NLO calculation of the n— jet 
cross section was used as input), while all higher orders in should be given by the parton 
shower rcsummation (possibly with improvements). This problem can be overcome by 

(1) Using tree-level matrix elements only as seeds for higher-order corrections, i.e. includ- 
ing the full resummation in tree-level events, and safeguarding that the weighting of 
n-jet tree-level configurations does not introduce 0(as"^^)-terms. 

(2) Defining an NLO cross section for n-parton states that does not include n+1 resolved 
jets, and making sure that no (uncontrolled) 0(as"'"^)-corrections are introduced by 
the NLO calculation. 

(3) Adding the corrected tree-level and the NLO events. 

This means that we have to decide how to define NLO cross section for n-parton states that 
do not include n + 1 resolved jets. We call an NLO calculation "exclusive" if it contains 
weights for n-jet phase space points that include Born, virtual and unresolved real correc- 
tions, where the resolution criterion is defined by exactly the same function as the merging 
scale. If all real emission corrections are projected onto n-jet kinematics, the calculation 
will be called "inclusive" instead. It is feasible to make an inclusive calculation exclusive by 
introducing explicit counter-events that arc distributed according to the resolved-emission 
contribution, and subtracting these events from events generated according to the inclusive 
NLO cross section. 

The points (1) - (3) will schematically lead to an algorithm of the form 

• Reweight n-resolved-jet tree-level events with weight used in the tree-level merging 
scheme. 

• Subtract the ©(a") and 0{ag~^^) terms introduced by this prescription from the 
tree-level events. 

• Add NLO-weighted n-resolved-jet events. 

Many variations of this basic form are possible. The first conceptual paper on NLO merging 
[29] for example advocated subtracting C'(a^"'"^)-terms from the NLO cross section. This 
is also the case in the MINLO NLO matching scheme [30], and the NLO merging scheme 
introduced for aMC@NLO [28]. MEPS@NLO [26,27] exerts full control over the NLO 
calculation to avoid point (1). We hope that the near future will bring detailed comparisons 
of all these schemes. 

Moving beyond tree-level merging prescriptions has long been regarded the next crucial 
step in background simulations for the LHC. The aim of this article is to present a com- 
prehensive guide to NLO merging schemes in PYTHIA8 [31]. We will present two different 
NLO merging schemes, choosing to generalise both the CKKW-L and UMEPS methods. 
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We will refer to the NLO version of CKKW-L as NL^, since this is an extension of the NLO 
merging scheme for e'''e~-coIlision presented in [25] to hadronic collisions, now also allow- 
ing for POWHEG input. The virtue of this method is its relative simplicity in combining 
NLO accuracy with PS resummation. However, it inherits violations of the inclusive cross 
section from CKKW-L. Cross section changes are the result of adding higher-multiplicity 
matrix elements ~ containing logarithmic terms that are beyond the accuracy of the parton 
shower, and can thus not be properly cancelled. At tree-level, this issue was resolved by the 
UMEPS method. Thus, we believe that a NLO generalisation of UMEPS, which also can- 
cels new logarithmic contributions appearing at NLO, is highly desirable. This method will 
be coined UNLOPS for unitary next-to-leading-order matrix element and parton shower 
merging^ . 

This publication is split into a main text and a large appendix section. The main 
text should be regarded as an introduction of the methods, while all technical details and 
derivations are collected in appendices. Sections in the main body are intended to give 
an overview of NL^ and UNLOPS, and explain some benefits with simple examples. The 
appendices are aimed at completeness, and in principle allow an expert reader to implement 
our methods in detail. As such, the appendices can also be considered a technical manual 
of the Pythias implementation. 

We begin by reviewing the CKWW-L method (section |2.1| ) and the UMEPS improve- 
ment (section |2.2| ). Then, we move to NLO merging methods (section |3|), discuss NL^ in 



section |3.1|, and describe UNLOPS in section 3.2. After these, we show the feasibility of 



the NLO merging schemes by presenting results for W-boson production (sections [4.1| and 



and for H-boson production in gluon fusion (section O). Then, we end the main text 
by concluding in section |5[ 

In appendix ^ we discuss some of the prerequisites that we need in order to derive our 
merging schemes, such as our choice of merging scale ( [A.l| ), the form of NLO input events 
that is required for NLO merging in PYTHIAS ( |A.2| ), a detailed description of the notation 
we use ( [A.3| ) and an outline of how the POWHEG-BOX program can be used to produce 
the desired NLO input ([A^). 

All technical details on how weights and subtraction terms are generated is deferred 
to appendix]^, ending in a summary (appendix |B.4D . From there, we move to a motivated 
derivation of the general NL^ method in appendix which can also be understood as a 
validity proof. The corresponding derivation of UNLOPS is given in appendix to which 
a comment on pushing this method to NNLO is attached (section p.l| ). We finally discuss 
the addition of multiparton interactions to NLO-merged results in appendix |^. 

Before moving to the main text, we would like to apologise for the inherent complexity 
of NLO merging methods. Also, we would like to affirm that in PYTHIAS, intricacies 
are handled internally, so that with reasonable input, producing NLO-merged predictions 
should not be difficult. The schemes described in this paper will be part of the next major 
Pythias release. 



^While finishing this article, a conceptual publication [32] was presented, which discusses similar meth- 
ods. 
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2. Tree- level multi-jet merging 

Since the methods presented in this article are heavily indebted to tree-level merging meth- 
ods, we would like to start with a brief discussion of CKKW-inspired schemes. Let us first 
introduce some technical jargon. 

We think of the results of fixed-order calculations as matrix element (ME) weights, 
integrated over the allowed phase space of outgoing particles. In the following, we will 
refer to a phase space point (i.e. a set of momentum-, flavour- and colour values for a 
configuration with n additional partons) as a n— parton event, n— parton state or simply 
S^n- Let us use the term "n-resolved-jet phase space point" (often shortened to n-jct point, 
n-jet configuration or n-jet state) for a point in the integration region for which all n partons 
have a jet separation larger than a cut tms- We further classify any configuration of hadronic 
jets by the number of resolved jets from which it emerged. The goal of merging schemes 
is to describe configurations with n hadronic jets with n-jet fixed-order matrix elements, 
meaning that the distribution of hadronic jets is governed by the "seed" partons in the 
n— jet phase space point, while parton showers only dress the seed partons in soft /collinear 
radiation. We will often use the notation O (S^nj) to indicate that the observable O has 
been evaluated on configurations containing n resolved jets. 

Tree-level merging schemes have to ensure that an n-resolved-jet configuration never 
evolves into a state with n -I- 1 well-separated hadronic jets. This can be achieved by 
applying Sudakov factors to the ME input events and by vetoing emissions above Ims- To 
capture the full parton shower resummation, it is common to also reweight the input events 
with a running coupling. Because no n-resolved-jet event evolves to an (n + l)-rcsolvcd- 
jet state, it is possible to add all contributions, and combine the tree-level description of 
well-separated jets with the resummation of parton showers, which can then be processed 
by hadronisation models. 

To summarise, tree-level merging is realised by calculating tree-level- weighted n-jet 
phase space points for up to N additional jets, reweighting these events, guaranteeing that 
ME events do not fill overlapping regions of phase space, and combining the different event 
samples for predictions of observables. It is not reasonable to limit observable predictions 
to only include configurations with up to N additional resolved jets. Instead, the parton 
shower is used to generate resolved jets for all multiplicities n > N, fov which no tree- 
level calculation is available. This is accomplished by not restricting the emissions off the 
highest multiplicity (iV-jet) ME state to unresolved partons only. 

There are in principle different ways how to combine the reweighted samples in tree- 
level merging. CKKW-inspired methods use an additive scheme, while unitarised MEPS 
opts for an add-subtract prescription. In the following, we will briefly discuss CKKW-L 
and UMEPS. 

2.1 CKKW-L 

The CKKW-L method [2,3,5] imposes tree-level accuracy on the parton shower description 
of phase space regions with n < N well-separated partons. To this purpose, tree-level- 
weighted phase space point are generated in the form of Les Houches Event files [33] . The 
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cross section of producing a state 5+„ with n partons has to be regularised by a cut ^ms 
on the momenta of the partons. In CKKW-L, any (collection of) cuts that regularise the 
calculation are allowed, tms is commonly called merging scale. 

The events with n = 1, . . . ,N additional partons will then be reweighted to incorporate 
parton shower resummation. The parton shower off states S'+„<7v is further forbidden to 
generate radiation that passes the cut ^ms- Reweighting with no-emission probabilities, 
and ensuring that parton shower emissions do not fill phase space regions for which events 
are available from other ME multiplicities, will guarantee that there is no double-counting 
between events with different number of partons in the ME calculation. In this publication, 
we will use the minimal parton shower evolution variable p as regularisation cut, and hence 



denote the merging scale by Pms- This choice is discussed in appendix [A.l . 

The full CKKW-L weight to make n-parton events exclusive, and minimise the depen- 
dence on Pms, is given by 



i=l 



X 



nSSj ^ \Jl^S+,^Axi-l,Pi-l,Pi)j X '^S+AXn, Pn, Pms) (2-1) 



n 

X 

= 1 



(^s{pi)x+_J+_-^^{x+_^,Pi.i)x-_J-_-^ix-_-^,Pi.l) , , 

-^l■S+^^l[X^-l, Pi-l, Pi) 



n as{pR) ^t-ift-li^t-i'Pi) ^i-^_fi-i{^i-i,Pi) ^^'^+* 
X '^S+AXn, Pn, Pms) , (2-2) 



where pi are reconstructed emission scales. The first PDF ratio in eq. (|2.1| ) means that the 
total cross section is given by the lowest order Born- level matrix element, which is what 
non- merged PYTHIA8 uses. The PDF ratio in brackets comes from of the fact that shower 
splitting probabilities are products of splitting kernels and PDF factors. The running of 
as is correctly included by the second bracket. Finally, the event is made exclusive by 
multiplying no-emission probabilities. The PYTHIA8 implementation reorders the PDF 
ratios according to eq. (^), so that only PDFs of fixed flavour and x- values are divided, 
thus making the weight piecewise numerically more stable. This will also later be useful 
when expanding the CKKW-L weight. For the highest multiplicity, the last no-emission 
probability li^s+Ni^N, Pn, Pms) is absent to not suppress well-separated emissions for which 
no ME calculation is available. 

The calculation of the CKKW-L weight is made possible by using a parton shower 
history. Parton shower histories are crucial for all merging methods, so it is necessary to 
elaborate. The matrix element state, S^n, (read from a LHE file) is interpreted as the 
result of a sequence of PS splittings, evolving from a zero-jet state, S+o, to a one-jet state, 
5+1, etc. until the the state S+n-i splits to produce the input S+n- All splittings occur at 
associated scales pi,. . . ,pn- A parton shower history (short PS history) for an input state 
S^n, i-e. a sequence of states, S+q, . . . , S+n, and scales pi, . . . , pn, is constructed from the 
input event by inverting the parton shower phase space on S+n. This means that in a first 
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step, we identify partons in S'+„ that could have resulted from a splitting, recombine their 
momenta, flavours and colours, and iterate this procedure on . . . , S+o- Clearly, 

there can be many ways of constructing such a history "path". Indeed, we construct all 
possible parton shower histories for each input state S+n, and then choose one history path 
probabilistically, using the product of PS branching probabilities as discriminant. More on 
this matter can be found in [5]. 

A CKKW-L merged prediction for and observable O is obtained by adding all contri- 
butions for fixed numbers of resolved jets 0{S+nj), given by the reweighted (i.e. exclusive) 
the S-^-n events, for all multiplicities n = 1,...,N. Using the symbol B„ for the fully 
differential n-parton tree-level hadronic cross section, the prediction for O is given by 



where we have used the symbol S+nj to indicate states with n resolved jets, resolved 
meaning above the cut /Jms as defined by the merging scale definition. The contribution of 
states with more than N resolved jets is included by allowing the parton shower to produce 
emissions above the merging scale when showering the A^-jet ME events. 

2.2 UMEPS 

The idea of unitarised matrix element + parton shower (UMEPS) merging [6] is to sup- 
plement CKKW-L merging with approximate higher orders for low-multiplicity states, in 
order to exactly preserve the n-jet inclusive cross sections. In UMEPS, events with addi- 
tional jets, which are simply added in CKKW-L, are also subtracted, albeit from lower- 
multiplicity states. This subtraction is motivated by the mechanism for how non-corrected 
parton showers would preserve the inclusive cross section. The contribution for a jet being 
emitted off S-^q at scale p, for example, is cancelled with contributions for no jet being 
emitted between Pmax and p. UMEPS makes this cancellation explicit by constructing 
subtraction terms through integration over the phase space of the last emitted jet. The 
guiding principle is "subtract what you add" : If n— parton events are added, those events 
should, in an integrated form, be subtracted from (n— 1)— parton states. Improvements in 
multi-jet observables are retained, since integrated n-parton events and "standard" events 
contribute to different jet multiplicities. 

In UMEPS, Les Houches events with (initially) n— partons are reweighted with 

fn i^n -iPn) fn (^n :Pn) 

xt (xt ,Pf) X~ f~ {Xn ,fJ-F) 




(2.3) 



(2.4) 



^S+i-i{Xi-l, Pi-1, Pi) 



(2.5) 
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This is the CKKW-L weight Wn without the last no-emission probabihty Ils+„ {xn, Pn, Pms) 
(i.e. for the highest multiphcity N: w'j^ = wn). As before, we make use of a PS history 
to calculate this weight. Denoting the tree- level differential n-parton cross section by B„, 
and introducing the notation 

Bnw'^ = B„, and ^ d""™./) B„u;; = j^^n^m , (2.6) 

we can write the UMEPS n-jet merged prediction for an observable O as 



(O) = I #o|0(S+o,) 



Bo — / Bi^o 

s 



B 



2-j-O 



B 



Bi 



B 



2-)-l 



B 



+ 



+ 
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J s 



N 
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E 

n=0 



0{S, 



Bn 



N 

E 

;=n+l 



(2.7) 



Many parts of standard CKKW-L implementations can be recycled to construct UMEPS 
predictions. The letter s on the integrals in the samples Bn+i^n indicates that the 
integrated states can directly be read off from intermediate states in the parton shower 
history. If a one-particle integration includes revoking the effect of recoils, it is possible 
that the state after performing one integration contains unresolved partons. In this case, we 
decide to perform further integrations (as indicated by the integration measure in 2.6), until 
the reconstructed lower-multiplicity state involves only resolved jets. Multiple integrations 
will include the effect of /OMs-unordered emissions into the description of lower-multiplicity 
states. We think of these (pMs-unordered, sub-leading) contributions as improvements to 
a strictly ordered parton shower. 

It might however not always be possible to find any parton shower histories that will 
permit at least one integration. If the flavour and colour configurations of a -l-n-parton 
phase space point cannot be projected onto an "underlying Born" configuration with n — 1 
partons, we call the parton shower history of the phase space point incomplete [5]. The 
existence of configurations with incomplete histories is reminiscent of which particles are 
consider radiative partons, meaning that if W-radiation were allowed, a history 



cc 



udW" 



cc 



uu 



is possible, while otherwise, the history of cc — >■ udW is incomplete. Note that the effect 
of incomplete configurations on the cross section is minor, since such contributions are 
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related to flavour changes of fermion lines through radiation. Configurations with incom- 
plete histories are not regarded as corrections to the lowest-multiplicity process, and will 
be treated as completely new process. Therefore, we will not (and cannot) subtract config- 
urations with incomplete histories from lower- multiplicity states, which leads to marginal 
changes in the inclusive cross section. 

2.3 Getting ready for NLO merging 

Following appendix |A.1| , we define the merging scale in terms of the shower evolution 
variable, thus putting tus = Pms- We further rescale the weights Wn and w'^ by a K- 
factor, K = J B/ J B, to arrive at a better normalisation of the total cross section. This 
means the introduction of additional 0(03(^/2))— terms, which have to be removed later 
on. Appendix |b| discusses the generation of these if -factors, which were introduced in [25] 
to avoid dicontinuities across the merging scales. Note that we include if-factors only 
because we do not see a formal reason against rescaling. In this publication, we attempt 
to provide a general definition of our new NLO merging schemes, and thus include these 
factors. 

Parton showers make as a tunable parameter, so that e.g. as(Mz) is chosen to fit 
data as closely as possible. This means asps(Mz) used in the parton shower might not 
be the value asME(Mz) used in the matrix element calculation. We can recover a uniform 
ag-definition by shifting 

cuspsip) = asME{bip), (2.8) 

where bi might be take different values, 6/ or if asps(Mz) is different for initial and final 
state splittings. If asME(^//F/o) would then be used instead of as everywhere, a uniform 
ag definition would be recovered. For this paper, we choose Qsps(Mz) = asME(Mz) = 
as^PDp(Mz), i.e. fix the value of as(Mz) to the one used in the parton distributions. In the 
future, when developing a NLO tune, we will interpret asps(Mz) as a tuning parameter, 
so that we can check the influence of NLO merging on the (rather high) parton shower 
as value. For the results in this publication, we will drop the index "ME" on as, and 
understand hi = 1. Our starting point for NLO merging are n— parton samples reweighted 

by 

„,J jy' fn i^n tpn) fn iS^n iPn) 

— ^ ■ —FTTTZT 



n 

X 

i=l 

Wn = w'^Tis+A^n, Pn, Pms) , (2-10) 



n 



xtjtli-t.P^) --_,/-,(x-_,,p.) ^^S+.-AX^-l,P^-l,P^) 



(2.9) 



in the case of UMEPS and CKKW-L, respectively. When referring to the weight in UMEPS 
and CKKW-L we will from now on always allude to the weights including a ii'-factor. 

Since we aim at interfacing two different program codes - NLO matrix element gen- 
erators and parton shower event generators - we need to make sure that the output of 
one stage (i.e. the NLO ME generator) is completely understood, before using it as in- 
put for the event generation step. Thus, we require that all fixed-order calculations are 
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performed with fixed factorisation and renormalisation scale, since dynamic scale choices 
in the fixed-order calculation result in subtle changes in higher orders^. All higher-order 
terms due to ag— running and PDF evolution will be carefully taken into account in the 
merging algorithm by reweig hting with eq. (^) (or eq. (p^). 



3. Next-to-leading order multi-jet merging 

Before sketching the NLO merging schemes we want to present here, we apologise that the 
discussion is (even after shifting most technical details into appendices) unfortunately very 
notation- heavy. 

Multi-jet merging schemes act on exclusive fixed-order input. This, for example, means 
that all phase space points that are allowed in the evaluation of tree-level matrix elements 
with n outgoing partons correspond to configurations with exactly n resolved jets, and no 
unresolved jets. The resolution criterion is given by the minimal separation of jets, with the 
relative transverse momentum used as shower evolution variable defining the separation^. 

The idea of using exclusive inputs is adopted for NLO merging, however, the notion 
of exclusive cross section needs to be refined for next-to- leading order calculations: We 
consider an n— jet NLO calculation exclusive if the output consists of n— parton phase 
space points with weights that correspond to the sum of Born, virtual and unresolved real 
radiation terms, where by unresolved real emission, we mean that the additional emission 
does not produce an additional resolved jet. It is possible to amend the NLO merging 
scheme if the requirement that all real emission terms are unresolved is not met (see 



discussion about exclusive vs. inclusive NLO calculations in appendix A. 2 for details). 

For an NLO merging scheme it is however crucial that virtual and unresolved real 
contributions contribute to the same phase space points, since otherwise, it is not possible 
to guarantee an implementation that is independent of the infrared regularisation in the 
NLO calculation. This problem is solved in POWHEG and MC@NLO, where real-emission 
contributions are projected onto n— jet phase space points by integrating over the radiative 
phase space. In this article, we use the POWHEG-BOX program [19] as NLO matrix element 
generator^ . 

Note that we do not require any change in the NLO matrix element generator. It is 
acceptable to produce LHE output with only minimal cuts. The merging scale jet sepa- 
ration will then be enforced internally in PYTHIA8, meaning that after reading the input 
momentum configuration from LHE file, any event not passing the cut will be dismissed. 
Pythias itself can decide if the required number of resolved jets are found, thus rendering 
the input exclusive. 

The aim of this section is to briefly describe two NLO merging algorithms. Each 
description will be split into a more formal part, and an algorithmic section, with the goal 
of presenting an overview of the NLO merging prescriptions coined NL^ and UNLOPS. So 



The preparation of output of the POWHEG-BOX program [19] is outlined in appendix A. 4 



''See appendix 
^See appendix 



A.l for details. 



A.4 for details. 
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that the flow of the narrative is not overly cluttered with technicalities, we have shifted all 
details into appendices. We however wish to introduce the reader to the symbols^ 



n+l|w 



B 

^n+l\n- 
-^n+l|n' 

B„: 

Is Bn— 
Is Bn— >n 

Bn: 

Ig Bn— >n 



Tree-level matrix element for n outgoing partons. 

Sum of tree-level cross sections with n outgoing partons in the input ME 
events, after integration over the phase space of n — m partons. 
Sum of tree-level configurations with n + 1 partons with a definite corre- 
spondence to a n-parton tree-level matrix element. 
Virtual correction matrix element for n outgoing partons. 
Sum of infrared regularisation terms for n resolved and one unresolved 
parton. 

Sum of integrated infrared regularisation terms for n resolved and one un- 
resolved parton. 

Inclusive NLO matrix element for n outgoing partons, i.e. sum of Born, 
virtual and all real contributions as weight of n— parton phase space points. 
Exclusive NLO matrix element for n outgoing partons, i.e. sum of Born, 
virtual and unresolved real contributions as weight of n— parton phase space 
points. 

Inclusive NLO cross sections with n outgoing partons in the input ME 
events, after integration over the phase space of n — m partons. 
Exclusive NLO cross sections with n outgoing partons in the input ME 
events, after integration over the phase space of n — m partons. 
UMEPS-processed n-resolved-jet tree-level events. 

UMEPS-processed tree-level cross sections with initially n resolved jets in 
the input ME events, after integration over the phase space of n—m partons. 
Contribution A, with terms of powers and Qg removed. 
Contribution A, with only terms of power ctg and ctg retained. 



Appendix is intended to give more thorough explanations of the notation. Particularly 
the last two short-hands are helpful when isolating orders in Og. For example, we have 

= 



-2 



B2 

Bo 



-0,1 



V„ -I- In+l|n + / d^ra.d B„+i|„6 (/9ms " ^{S+n+l, p)) - Dn+l| 



Boi 



Wo 



Bn 



115+0(^0, /50,Pms) 



' PMS 



27r 



aE{outgoing} j 



at \incomingj J 



See appendix A. 2 for details 
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All details on the expansion of the tree-level weights can be found in appendix ^ We are 
now equipped for extending CKKW-L and UMEPS tree-level merging to next-to-leading 
order accuracy. 

3.1 NL^: CKKW-L at next-to-leading order 

The NL'^ prescription [25] in principle starts from the CKKW-L-weighted tree-level cross 
sections BnWn, adds events weighted according to the exclusive NLO cross sections B„, 
and removes approximate 0{a^{fj,ji)) and 0{al{fj,ji)) terms in the CKKW-L weight Wn- 
Since exclusive NLO samples are rarely accessible, we instead use the inclusive NLO cross 
section and generate explicit subtraction events by using higher- multiplicity tree-level 
matrix elements. For details of this choice, we refer to appendix |A.4| . 

All details about the derivation of the NL^ method can be found in appendix |C|. Here, 
let us assume the construction of NLO accuracy -|- parton shower higher orders is possible 
for configurations with exactly m resolved jets, and that the desired accuracy is achieved for 
any number of resolved jets m € {0, . . . , M}. On top of these NLO-correct multiplicities, 
NL^ allows the inclusion of tree- level matrix elements with n S {M -|- 1, . . . , A^} additional 
partons. The highest-multiplicity tree-level sample further allows the generation of more 
than N resolved jets, by allowing parton shower emissions to produce resolved partons. 
The complete result is then obtained by simply adding the partial results for each jet 
multiplicity. This means that the NL^ result for an observable O, when merging N tree- 
level, and M < N next-to-leading order calculations, is 

M 



m=0 

N 



-m,m+l 



B 



(3.1) 

n=M+l '' J J 

where the crucial change from CKKW-L (c.f. eq. (|2.4D ) is in the first line, where we add the 
exclusive NLO events and remove the corresponding ag-terms from the CKKW-L weight. 
From a technical point of view, it is often convenient to think of this in terms of processing 
the samples 

J for m<M (3.2) 

for m< M (3.3) 
for m < M (3.4) 
for M < n < A (3.5) 















— m,m+l 






Bm 




Sm 




- / Bm 


+ 1— 






J s 




Tn 




BnWn 





and writing simply 



(0) = E / 

m=0 

N 

+ E 



OiS- 



+mj) 



n=M+l 



(3.6) 
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The prediction for a generic observable can be obtained by calculating the result Oh{S^k)j 
measured for /c— jet phase space points, filling the histogram bin Ob with weight (or 
T'^/Vfc/Sfc, depending on the sample), and summing over all multiplicities k. 

The construction of the necessary weights is done with the help of a parton shower 
history and is detailed in appendix ^ Once the weights are calculated, further parton 
showering is attached. The shower off inclusive NLO events and phase space subtractions 
is started at the last reconstructed scale, and all emissions above Pms vetoed. This means 
that all higher-order terms above the merging scale are taken solely from the reweighted 
tree- level matrix elements, thus ensuring that the prescription preserves the parton shower 
description corrections beyond the reach of the NLO calculation. All samples have to be 
added to produce NLO-accurate m = 0, . . . ,M jet observables, with higher-order correc- 
tions given by CKKW-L. Details on how the weights of different samples are motivated, 
as well as a proof of NLO -|- PS correctness, are given appendix 

Here, let us illustrate how NLO accuracy is achieved for one particular jet multiplicity. 
For this, we examine the samples contributing to M— jet observables (where M is the 
highest multiplicity for which an NLO calculation is available). We start by analysing the 
0{a^) and 0(0^+^) contributions. We find 



M 



+ 



{O) 



M 



M+l 



0(5+m,)Vm 



M 



+ 



0{S+Mj)\BM + VAf + lM+l|Af + j d^rad (B 



0{S+mj)[Vm + Sm] 

M+1\M — ^M+1\m) 



M+l 

Ba/+i^a/ 



} 

(3.7) 



Thus, the description of M— jet states is NLO-correct. For M + 1— jet events, we have 



{0)m- 



M+l 



0{S-^Mj) Ba/ 



+1 



(3.8) 



providing tree-level accuracy. Both these facts mirror the NLO description of observables. 
Keeping only the next-higher powers 0{a^~^'^) above Pms, we see that 



{O) 



M 



+ 



{O) 

M+2 L 
0{S^Mj) BAf 



A/+1 



A/+2 



+ 



M+2 



M+2 



(3.9) 



WM 



+ 0(5'+A/+ii) Baz+i 



WM+1 



+ 0{S+M+2j) Ba/+2 



For M — jet observables, only the reweighted M— parton LO matrix element contributes, 
while M + 1— jet observables are described by reweighted M + 1— parton tree-level states. 
M-|-2— jet observables are determined by the M-l-2— parton tree-level prediction. These are 
the results of default CKKW-L. Thus, the method is NLO accurate for M— jet observables, 
and also retains exactly the resummation of CKKW-L in higher orders for M— and M + 
1— jet observables. 

3.1.1 NL^ step-by-step 

In the NL^ algorithm, we have to handle three classes of event samples: 
A: Inclusive next-to-leading order samples Vm for m < M resolved jets. 
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B: Tree-level samples for m < M resolved jets, and tree-level samples T„ for M < 
n < N jets. 

C: Tree-level samples §rn with initially m+ 1 partons, after integration over the radiative 
phase space of the (m -I- l)'th parton (for rn < M). 

Samples of class A are produced with the POWHEG-BOX program, by setting the minimal 
scale for producing radiation to Ecu- For calculations that need to be regularised, we use 
minimal cuts in POWHEG-BOX, and reject events without exactly the number of required 
jets internally in PYTHIA8. The samples of class A are processed in the most simple 
manner: 

A. I Pick a jet multiplicity, m, and a state S'+„t) according to the cross sections given by 

the (NLO) matrix element generator. Reject any state with unresolved jets. 
A. II Find all parton shower histories for S-^-o, ■ ■ ■ , S+rm and pick a parton shower history 
probabilistically. 

A. Ill Do not perform any rcweighting on S+m- 

A. IV Start the shower off S+m at the latest reconstructed scale pm- Veto shower emissions 
resulting in an additional resolved jet. Pms- 

A. V Start again from A.I. 

To amend that we have used inclusive NLO cross sections where we should have used exclu- 
sive calculations, we have to introduce samples of class C. The first step in the construction 
of these samples is to generate tree-level weighted events with 1 < m < M + 1 partons 
above Pms- Then, 

C.I Pick a jet multiplicity, m + 1, and a state S+m+i, according to the cross sections 
given by the (LO) matrix element generator. Reject any state with unresolved jets. 

C.II Find all parton shower histories for S+q, . . . , S-^-m, S+m+i, and pick a parton shower 
history probabilistically. Replace (S+^+i with the -S+m of by the chosen history^. 

cm Weight S+m with -1. 

CIV Start the shower off S+m at the latest reconstructed scale pm- Veto shower emissions 
resulting in an additional resolved jet. 

CV Start again from C.I. 

Higher orders in Os (in the CKKW-L scheme) are introduced by including events of class 

B. Again, tree-level weighted events for < n < A'^ partons are needed as input. Then, 

B. I Pick a jet multiplicity, n, and a state S-^-n, according to the cross sections given by 

the matrix clement generator. Reject any state with unresolved jets. 

B.II Find all parton shower histories for S+q, . . . , and pick a parton shower history 
probabilistically. 

B.III Perform reweighting: 

'^We do not apply any further action if S+m contains unresolved jets in NL^, in contrast UMEPS (or 
UNLOPS). 
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B.III.l If n > M, weight with Wn, as would be the case in CKKW-L. 
B.III.2 If n < M, weight with {wn — [wuIq — [wnli}- 
B.IV Start the shower off 5'+n at the latest reconstructed scale pn- 
B.IV.l If n = A^, allow any shower emission. 

B.IV. 2 If n < A^, veto shower emissions resulting in an additional resolved jet. 
B.V Start again from B.I. 

All samples of all classes are finally added to produce the M-NLO-jet- and A^-LO-jet- 
merged prediction. Both the samples B and C require tree-level input, i.e. the input events 
for C-samples can be also be used as input for B-samples. In total, the PYTHIA8 imple- 
mentation requires M NLO-weighted Les Houches event files, and tree-level-weighted 
files as input, but some of the tree-level input files need to be processed twice. 

Due to the ubiquity of multiparton interactions (MPI) in hadronic collisions, we are 
still far from a full event description at the LHC, even after combining multi-jet calculations 
and parton showers. How MPI can be attached to NL^ is discussed in appendix]^. 

3.2 UNLOPS: UMEPS at next-to-leading order 

Although NL^ accomplishes a merging of multiple NLO calculations to the specified ac- 
curacy, it inherits the merging scale dependence of the inclusive lowest multiplicity cross 
section from CKKW-L. For lack of a better term, we will refer to changes in the inclusive 
cross section as "unitarity violations". When including additional jets in W-boson pro- 
duction, unitarity violations enter at the same order in as as e.g. the NLO corrections to 
W -|- j— production. Even if changes of the inclusive cross section are generally small as 
long as the merging scale is not set too small, it is not clear how much of the shape changes 
we observe are really due to not cancelling logarithms. Thus, we want to promote UMEPS, 
where these unitarity violations are absent [6], to NLO accuracy as well. 

Extending UMEPS to include multiple NLO calculations is slightly more involved 
than the CKKW-L case. The complete method is derived in appendix |D]. In a sense, NL^ 
and UNLOPS are complementary: NL^ is, in the accuracy claimed by the method, easily 
applicable to any number of jets, while UNLOPS aims at higher accuracy for the dominant 
low multiplicities^. The strategy to extend UMEPS to NLO accuracy is similar to NL^. 
We remove any approximate O (ag (/ij?)) and O [al (/ii?)) terms in the UMEPS weighting 
procedure, and simply add the correct NLO result. To disturb the description of higher 
order contributions as little as possible, we only cancel those terms of the UMEPS weight 
that would have a better description in the NLO matrix element. 

The UNLOPS method aims to move beyond UMEPS not only in terms of fixed-order 
accuracy for multiple exclusive n— jet observables, but also in the description of higher 
orders in low-multiplicity states. This is a direct consequence of requiring unitarity, i.e. 
that the inclusive cross section be fixed to the zero-jet NLO result. In the spirit of UMEPS, 

*In fact, the UNLOPS zero- and one-jet NLO merging presented here can easily be promoted to a NNLO 
matching scheme, as outlined in appendix O 
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this means that once we want to add a one-jet NLO calculation, we have to subtract its 
integrated version from zero-jet events. Similarly we need to remove the O (a° (hr)) and 
O ifJ-R)) in the UMEPS tree-level weights, not only for the one-jet events but also for 
the corresponding subtracted zero-jet events. In this way we ensure that the inclusive 
zero-jet cross section is still given by the NLO calculation and we will also improve the 
©(ag)— term of exclusive zero-jet observables. 

The UNLOPS prediction for an observable O, when simultaneously merging inclusive 
NLO calculations for m=0, . . . , M jets, and including up to N tree-level calculations, is 
given by 



M-l 



m=0 



{0)= / #0 / • • • / OiS+^j) \ B™ + 



B, 



M „ M 

J ^i^rn - ^ 
i=m+l ^ i=m+l 



-in,m+l 



B,; 



N 



- ^ r- 



N 



0{S+Mj){ Bm + 



i=M+l 

+ Yl [d^O [■■■ [ 0{S+nj) Yl I ^i^n | 

n=M+l J J J y i=n+l ) 



(3.10) 



Here we see, in the first line, the addition of the and the removal of the O (a^ (mk)) 
and O (a™^^ (M-r)) terms of the original UMEPS B^ contribution. On the second line we 
see the subtracted integrated B^+i term to make the m-parton NLO-calculation exclusive 
and the corresponding O (a™ {fJ-n))- and O (a™"^^ (/i/j)) -subtracted UMEPS term together 
with subtracted terms from higher multiplicities where intermediate states in the clustering 
were below the merging scale. The third line is the special case of the highest multiplicity 
corrected to NLO, and the last line is the standard UMEPS treatment of higher multiplic- 
ities. 

The full derivation of this master formula is given in appendix |D[ where we also discuss 
the case of exclusive NLO samples and explain the necessity for subtraction terms from 
higher multiplicities. We will limit ourselves to including only zero- and one-jet NLO 
calculations in the results section. For the sake of clarity we will thus only discuss this 
special case here. For this case, the UNLOPS prediction (when including only up to two 
tree-level jets) is 



(O) = / d<Ao|0(5+o,) I^Bo - 

+ // 0{S+2j)B2 



B 



1^0 



Bi 



-1,2 



Bi^o 

B2^1 



-1,2 Js 



B 



2^0 



(3.11) 
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In an implementation, this is conveniently arranged in terms of the samples 



Vn 



B2W2 
Bi 

Bn 



B 



-1,2 























.} 



Bo + Vo + Il|0 + / d$rad (Bi|o " Di|o) 



Vi = Bi = Bi + Vi+I 



211 



lo 
Lo 



/ Ba^i^ - /b 

J s J s 



d^rad (B 



211 



D 



211 



2^0'J^2 



Bi^o 
Bi^o 



-1,2 



Bi^o I w'l 



(3.12) 
(3.13) 

(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 



meaning that we have two tree-level samples (Bi, B2), two NLO samples (Vo, Vi), two 
subtractive samples {Iq, Ii) and one integrated NLO sample (Lo). The prediction is formed 
by reading tree- level input events (for Bi, B2, Iq and Ii), or inclusive NLO input (for Vq, 
Yi and Lq), generating the necessary merging weights, and filling histogram bins with the 
product of matrix element and merging weight. For technicalities on the generation of the 
weights, we refer to appendix ^ 

In the inclusive cross section, it can immediately be checked that all contributions 
except zero-jet NLO terms cancel exactly, meaning that the inclusive cross section is given 
by the zero-jet NLO cross section. As in UMEPS however, we have to accept marginal 
changes of the inclusive cross section in the presence of incomplete histories, i.e. when it is 
not possible to regard n-jet states as corrections to n — 1-jet states, because no underlying 
Born configuration exists, (see discussion at the end of section [2^ ) . The contribution from 
such configurations is, for the results presented in this publication, numerically insignificant. 

Let us turn to the UNLOPS description of exclusive observables. Only looking at 
zero-jet observables in eq. ( |3.11| ), we see 

{O)o = J #oO(S+o,) I B0 + V0 + I110 + J (i^'rad[Bi|oG(pMs-t(5+i,p))-Diio 



Bi 



^ Vi + I2I1 + j d^rad (B2II - D2I1) 

>0 I w'l - w'l - w[ \ - B2^0'«^2 
^ LJO LJIJ J g 



(3.19) 



where we have, between the first and second lines, cancelled the tree-level contribution 
of J"^ Bi^o with the resolved real-emission term J"^^^ (i$radBi|o appearing in Bo. When 
extracting only the O [a^ (/i/j)) and O [al (/i/j)) terms, this gives the contribution of the 
exclusive NLO matrix element, i.e. of tree-level, virtual correction and unresolved real 
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contributions. At O (fiR)), we have 

{O)o\^ = J d(PoO{S+oj)i^- V1 + I211 + y d$rad (B211 - D211) -^B2^o| (3-20) 

The first group of terms in the curly brackets gives an approximation of NNLO corrections, 
since in a NNLO calculation, all logarithmic terms in the NLO +1— jet calculation are 
removed by two-loop and double-real terms. Conversely, we should be able to include the 
correct logarithmic terms of two-loop and double unresolved terms by integrating over the 
jet in the +1— jet NLO calculation. The last term in curly brackets is sub-dominant and 
corresponds to emissions that are unordered in pMs ^■ 

Examining the O (ag {/j^r}) contributions, we are left with 



B 



1^-0 



B2^-0 



W2 



(3.21) 



This is simply the parton shower approximation, amended with a term corresponding 
emissions that are unordered in Pms- 

Let us move on to the discussion of one-jet observables. If we use the fact that we can 
cancel the contribution of two resolved real-emission jets in Bi by the O (oig {/j^r)) term in 
/^B2^i, we find 



{0)i ^ + 



(3.22) 



/r ( fPMs ^ 

dct>o j 0{S+ij) I Bi + Vi + I2I1 + j d$rad (B2II - D2I1) \ 



Thus, the method describes one-jet observables with NLO accuracy. The 
O (og {iir)) —term is given by 



Bi 

















(3.23) 



which is simply the UMEPS-improved parton shower approximation. In conclusion, we find 
the method is NLO-correct, improves the logarithmic behaviour of zero-jet observables, and 
otherwise includes the parton shower resummation of the UMEPS procedure. 

3.2.1 UNLOPS step-by-step 

As a complication on top of NL^, UNLOPS requires four classes of events. We will step- 
by-step formulate the method for including M inclusive next-to-leading order calculations, 
combined with N tree-level matrix elements. This means that we need to handle the 
samples 

A: Inclusive next-to-leading order samples for m resolved jets. 



^This term already appears in UMEPS. 



- 18 - 



B: Tree-level samples B„ for n < N resolved jets. (There is no zero-jet tree-level sample.) 

C: Tree-level samples I„ with initially m partons, after integration over the radiative 
phase space of the one or more emissions, as required by the UMEPS method. 

D: Next-to-leading order samples with initially m resolved jets, after integration over 
the radiative phase space of the emission. 

Samples of class A are produced with the POWHEG-BOX program, exactly as in NL^. The 
POWHEG-BOX output files are then processed: 

A.I Pick a jet multiplicity m, and a state S+m, according to the cross sections given by 
the (NLO) matrix element generator. Reject any state with unresolved jets. 

A.II Find all parton shower histories for iS+o, • • • , 5'+^, and pick a parton shower history 
probabilistically. 

A.III Do not perform any reweighting on 5+^. 

A.IV Start the shower off S+m at the latest reconstructed scale pm- Veto shower emissions 

resulting in an additional resolved jet. 
A.V Start again from A.I. 

This is exactly the treatment we already know from NL^. To ensure that the inclu- 
sive (lowest-multiplicity) cross section is not changed, we need to subtract the integrated 
variants I^m+i of the (m + l)-jet NLO calculation, i.e. introduce the samples of class D. 
This will also remedy the fact that we have used the inclusive m-jet NLO cross sections, 
while we should have used exclusive B^ input for Ym- As starting point, we use the 
Bm-i-i— distributed event sample (i.e. the same input as for V,„+i). Then 

D.I Reject any state with unresolved jets. 

D.II Find all parton shower histories for 5+0, . . . , 5+^+1) and pick a parton shower history 
probabilistically. Replace S+rn+i with the 5'+^, or the first state S+i with all I < 
m partons above the merging scale, (lower multiplicity states are taken from the 
intermediate states of the chosen PS history) 

D.III Weight S+rn with -1. 

D.IV Start the shower off 5+m at the latest reconstructed scale pm- Veto shower emissions 
resulting in an additional resolved jet. 

D.V Start again from D.I. 

To add the UMEPS resummation to these samples (and correct that we have used Bj^f 
events rather than exclusive Bm input for Ym), wc include samples of class C. These are 
generated from the (n -|- l)-jet tree-level samples, by following the steps 

C.I Pick a jet multiplicity, n + 1, and a state according to the cross sections given 

by the matrix element generator. Reject any state with unresolved jets. 

C.II Find all parton shower histories for S+q, . . . ,-S'+„, and pick a parton shower history 
probabilistically. 

cm Perform reweighting: 
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C.III.l If n + 1 > M, weight with as would be the case in UMEPS. 

C.III.2 If n + 1 < M, weight with - - J. 

CIV Replace S+n+i with the or the first state 5+/ with all / < n partons above the 
merging scale (lower multiplicity states are taken from the intermediate states of the 
chosen history). Start the shower off S+n at the latest reconstructed scale pn- Veto 
shower emissions resulting in an additional resolved jet. 

C.V Start again from C.I. 

The last contributions we have to include are reweighted tree-level samples, i.e. 
events of class B. There is no zero-jet tree-level contribution in UNLOPS, since the 
O (ag {fiR,)) —term is already included by Vq. Samples for class B are generated very 
similar to events of class C, with no "integration step" required for class B: 

B.I Pick a jet multiplicity, n > 0, and a state S+n, according to the cross sections given 
by the matrix element generator. Reject any state with unresolved jets. 

B.II Find all parton shower histories for S'+o, . . . , 5'+„, and pick a parton shower history 
probabilistically. 

B.III Perform reweighting: 

B.III.l If n > M, weight with would be the case in UMEPS. 

B.III.2 If n < TV, weight with {w[ - [w[]q - 

B.IV Start the shower off S'+„ at the latest reconstructed scale /)„. 

B.IV.l If n = A^, allow any shower emission. 

B.IV. 2 If n < A^, veto shower emissions resulting in an additional resolved jet. 

B.V Start again from B.I. 

Note that although the UNLOPS procedure is more complicated than NL'^, no addi- 
tional user input is required: PYTHIA8 only needs M inclusive NLO event samples and 
A^ — 1 tree-level event files, since A and B use the same NLO input, and the same tree-level 
input can be employed in both C and D. This concludes our discussion of NLO merging 
prescriptions. Information on how underlying event is added to our prescription is given 
in appendix 

3.3 Short comparison 

Before presenting results, let us pause and recapitulate the last section. We have presented 
two different NLO merging schemes, which differ in several ways 

NL3 UNLOPS 



• Generalisation of CKKW-L 

• Needs exclusive or inclusive NLO cal- 
culations as input. 



o Generalisation of UMEPS 
o Needs exclusive or inclusive NLO cal- 
culations as input. 
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• Straight-forward when moving to high 
jet multiphcities. 

• Changes the inclusive NLO cross sec- 
tions. 

• Reproduces the logarithmic behaviour 
of the PS in zero-jet observables. Does 
not fully cancel sub-leading logarith- 
mic enhancements of higher multiplicity 
NLO calculations. 

• Produces negative weights. 



o Less transparent when moving to high 
jet multiplicities. 

o Preserves the NLO inclusive cross sec- 
tions. 

o Explicitly cancels logarithmic enhance- 
ments, has improved logarithmic be- 
haviour in low- multiplicity jet observ- 
ables. 

o Produces even more negative weights. 



At this point, we will not make comparisons with other NLO merging methods, but hope 
to be able to contribute to a thorough comparison in a future publication. Here we will 
only make some brief remarks on the formal accuracy of our methods, compared to the 
ones presented in [26,27] (MEPS@NLO) and [28] (aMCONLO). AU of these methods rely 
on the introduction of a merging scale and it is relevant to investigate how the description 
of observables are affected by changes in this scale. In particular it is interesting to make 
sure that the NLO-correctness of the of the methods are not spoiled by large logarithms 
involving the merging scale, L = In Even if the dependence on the merging scale 

vanishes to the logarithmic approximation of the shower (normally at best NLL), the sub- 
leading logarithmic dependence may become as large as the NLO correction which we want 
include. 

To exemplify (following the arguments of Bauer et al. [34,35]) we look at the inclusive 
n-jet cross section, which in all methods have been corrected to reproduce the NLO cross 
section, so it is exact to ©(a") and 0{al!'~^^). But if we look at the 0(a""'"^)-term there will 
be dependencies on the merging scale, which we can symbolically expand out in logarithms 
as ^""'"^(L^ + + L"^ + . . .), Even for a NLL-correct parton shower where the both the 
Q,n+2^4 g^j^j^ a""^^L^ terms will cancel exactly, we will have dependencies of the order 
Q,ra+22^2_ 'pjjjg nieans we that have to choose the merging scale such that OgL'^ ^ 1, to be 
sure we do not spoil the effect of the ©(a" ^^)-correction of the NLO calculation. 

For the MEPS@NLO method it was shown that it at most has a dependence of order 
Q,n+22^3 ^ijig]^ is colour-suppressed, but certainly has a dependence of order a"+^L^. For 
the aMC@NLO method we do not know of any formal analysis of the logarithmic correct- 
ness, but it is difficult to see how it could have avoid dependencies of order a"+^L^. Also 
in our NL^ method, where the dependence is given by the precision of the shower, it can- 
not be claimed that the dependence of order a"+^L^ is absent, as it has not been proven 
that the PYTHIA8 shower is formally NLL-correct. However, for our UNLOPS method, we 
explicitly conserve the inclusive NLO cross section, and the merging scale dependence is 
cancelled almost completely through our "subtract everything that is added" strategy. We 
say almost cancelled, as this is clearly an observable-dependent statement. From eq. ( 3.1C| ) 



we see that in order for the addition of a higher order matrix element Bf^. to completely 
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cancel for an inclusive n-jet observable, we require (symbolically) 



fe-i 



/ #0 [■■■[ On{S+kj)Bk = I] / ^'^0 /• • • / On{S+^j)Bk^i, (3.24) 

i=n 

which is clearly never an exact cancellation. There is also an implicit merging scale depen- 
dence here as, whether or not e.g. is projected into B3^2 and measured with 0{S+2j) or 
into B^^i and measured with 0{S+ij), depends on the merging scale. However, for small 
enough merging scales, this should not matter for collinear- and infrared-safe observables, 
and we do not expect any large logarithms of the merging scale to appear. Also we note 
that there are some n-parton states that cannot be projected down to a lower multiplicity 



state using parton shower splittings {incomplete states in section p.2[) as described in [6], 
where we also found that such diagrams give numerically very small contributions. 

UNLOPS also shares features with the LoopSim method [36,37], in particular the use 
of an integrated version of one-jet NLO calculations. However, we cannot cancel logarithms 
of the form In (j^^^^ which arise by soft /collinear W-radiation, because we do not allow 
an integration over the (radiated) W-boson. The study of such "giant i^-factor effects" is 
postponed until a full electroweak shower becomes available in PYTHIA8. 

Finally it should be noted that NLO merging methods can be useful even if only 
the NLO calculation for the lowest multiplicity is available. Since an NLO merging scheme 
consistently splits the real emission into unresolved and resolved parts by defining a merging 
scale Pms) and uses the same definition to separate states with two resolved jets from states 
with one resolved and several unresolved jets, any NLO calculation can be improved by 
merging further tree- level calculations for additional jets. Such schemes go under the 
name of MENLOPS [38-40]. Promoting a NLO calculation to a MENLOPS prediction is 
straight-forward with our methods. 



4. Results 

The UNLOPS and NL'^ methods have been implemented in PYTHIAS, and will be included 
in the next major release version. In this section, we will present sample results for NLO 
merging with inclusive NLO calculations. The aim of this section is to affirm that the 
implementation in PYTHIAS is working smoothly. We do so by presenting results for W- 
boson production and Higgs (H) production in gluon fusion, when simultaneously merging 
zero and one additional jet at next-to-leading order with two additional jets on tree-level. 

All input matrix element configurations are taken from Les Houches Event Files. We 
use the following input: 

• W + 0, W + 1 and W + 2 at tree-level generated by MadGraph/MadEvent. 

• W + 0, W + 1 at NLO [41,42] generated by POWHEG-BOX (see appendix |AJ) . 

• H + 0, H + 1 and H + 2 at tree-level generated by MadGraph/MadEvent. 

• H + 0, H + 1 at NLO [43,44] generated by POWHEG-BOX (see appendix |AJ) . 
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• Fixed-order input was calculated with three values for fixed renormalisation scales 
and factorisation scales, 

- Central scales: ^lr = and = for W-production, = and 
HF = Ml = (125 GeV)2 for H-production. 

- Low scales: ^xr = (Mz/2)2 and ^xp = (Mw/2)2 for W-production, = (Mz/2)2 
and = (Mh/2)^ for H-production. 

- High scales: /ij? = {2M^f and ^ip = {2M^f for W-production, ijlr = {2M^f 
and = (2Mh)^ for H-production. 

In all Figures we will label curves generated from central scale input with cc, from 
low scale input with //, and from high scale input with hh. 

• CTEQ6M parton distributions and as{Ml) = 0.118. 

• The merging scale Pms is defined by the minimal PYTHIA8 evolution variable (see 
appendix |A.l| ). 

The value of a^i^M^) was set to match the Og-value obtained in the parton distributions 
used in the ME calculation. We use the same PDFs and Q;s(M2)-value in the parton 
shower evolution. For all internal analyses, we use fast jet-routines [45] to define jets. 
The momentum of the intermediate W-boson will, if required, be extracted directly from 
the Monte Carlo event. 

We will compare our results to the result of the POWHEG-BOX program for W+jet 
production. For these comparisons, we have generated default POWHEG-BOX output, fixing 
the renormalisation and factorisation scales, and regularising the Born configuration with a 
cut p^^parton = 5 GcV. To determine a shower starting scale for these POWHEG-BOX output 
events, we reconstruct all possible (including unordered) parton shower histories, choose 
one, and start the shower from the last reconstructed scale. No visible effects of using 
different options to choose history have been found. This is not the default interface to the 
POWHEG-BOX, which requires truncated showers if the scale definition on the POWHEG- 
BOX and the parton shower do not match. Appropriately vetoed showers are normally used 
instead in PYTHIAS, because no truncated showers are available. Since the scale definition 
in the POWHEG-BOX could change depending on the details of the implementation (being 
different for Catani-Seymour- and Frixione-Kunzt-Signer-based approaches), we do not use 
vetoed showers, and rather choose starting scale by constructing a parton shower history. 
For W-|-jet production, we found only insignificant differences between both methods. 

When taking ratios to default PYTHIAS (often given by the lower insets of figures), 

we rescale the PYTHIAS reference by K{^r,^f) = /BoI/ja'^f) ' This guarantees that we 
remove the variation of the normalisation of the inclusive cross section due to scale choices: 

jBo(Ma,MF)^^^^^° '^"'^^'^ ^ (g)NLO merged 

jB;;Q7^;^(C)pythia8 ^(/^i?, ^f) • (0)pythia8 

The variation of the overall normalisation will otherwise obscure interesting features. For 
Higgs production in gluon fusion for example we will compare merged curves generated with 



-23- 




20 40 60 80 100 120 140 160 ISO 20 40 60 80 100 120 140 160 180 

Pll [GeV] Pj^, [GeV] 



Figure 1: Transverse momentum of the hardest jet, for W-boson production in pp collisions at 
Eqm = 7000 GeV, when merging up to two additional partons at LO, and zero and one additional 
parton at NLO. Jets were defined with the /c^-algorithm, with k^^min = fO GeV. Multi-parton 
interactions and hadronisation were excluded. We choose to order the contributions by the jet 
multiplicity in the states on which further the showering is initiated. Left panel: Results of the 
NL'^ scheme. Right panel: Results of the UNLOPS scheme. 
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Figure 2: k±- separation dgi of the first jet and the beam, for W-boson production in pp collisions 
at EcM = 7000 GeV, when merging up to two additional partons at LO, and zero and one additional 
parton at NLO, for three different merging scales. Jets were defined with the fcj_-algorithm, by 
clustering to exactly one jet. Multi-parton interactions and hadronisation were excluded. Left 
panel: Results of the NL'^ scheme. Right panel: Results of the UNLOPS scheme. 

HR = 2Mz and fip = 2Mh (labelled hh) to PYTHIA8, multiplied by with i^(2Mz,2MH). 
For central scales, we would use i^(Mz,MH). 

4.1 W-boson production 

Let us start by discussing results for W-boson production, when combining inclusive NLO 
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Figure 3: Comparison of using exclusive and inclusive NLO input for W-boson production in pp 
collisions at i?cM = 7000 GeV, when merging up to two additional partons at LO, and zero and 
one additional parton at NLO. Curves labelled "inc" are produced with the B-prescription, while 
"exc" indicate a generation with B-input. The lower inset shows the deviation from PYTHIA8. 
The band labelled "POWHEG W+jet" is given by the envelope of varying the renormalisation 
scale in the POWHEG-BOX program between iMz, . . . , 2Mz, and the factorisation scale between 
iMw,...,2Mw. 

calculations for W + 0- and W + 1 parton with the PYTHIA8 event generator. This section 
is intended mainly for validation, and we will thus switch off multi-parton interactions 
and hadronisation. We present results for both NL^ and UNLOPS. Our preferred method 
is UNLOPS, since the inclusive cross sections are there handled more consistently. By 
showing the results for both NL^ and UNLOPS, we hope to convey a rudimentary idea of 
the effect of potentially problematic logarithmic enhancements in standard observables. 

Figure |l| shows the transverse momentum of the hardest jet in W-production at the 
LHC, for both NL^ and UNLOPS. The sum of the solid, coloured curves gives the full NLO 
merged result, i.e. the black line. The dashed curve is included only to illustrate that the 
hard tail of the p_Li-spectrum is dominated by the the W+jet NLO sample (labelled ^i), 
both in NL^ and UNLOPS, which is of course desired. This fact makes the pj_i-spectra of 
NL^ and UNLOPS very similar. 

Differences between NL^ and UNLOPS are expected in the intermediate- / low-scale 
regions. This is illustrated by Figure |2|, which shows the doi-distribution of the first jet^*^. 
Since UNLOPS explicitly preserves the W-production NLO cross section, the increase in 
the tail has to be compensated by decreasing contributions below the merging scale. The 
description at low doi in NL^ is, by construction, completely governed by the PYTHIA8 
result. 

Before continuing, we would again like to stress that we are using inclusive NLO cross 

^"The observable doi is very closely related to p±i , but avoids a fcx,min-cut in defining the jet, by clustering 
to exactly one jet. This allows to show the lowest scale features. 
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Figure 4: Transverse momentum of the hardest jet, for W-boson production in pp conisions at 
EcM = 7000 GeV, when merging up to two additional partons at LO, and zero and one additional 
parton at NLO. Jets were defined with the fc^-algorithm, with k±,rnin = 10 GeV. Multi-parton 
interactions and hadronisation were excluded. The lower inset shows the deviation from PYTHIA8. 
The band labelled "POWHEG W+jet" is given by the envelope of varying the renormalisation 
scale in the POWHEG-BOX program between ^Mz, . . . , 2Mz, and the factorisation scale between 
iMw, . • . , 2Mw. Left panel: Resuhs of NL^. Right panel: Resuhs of UNLOPS. 



section as input in this publication, as discussed in appendix A. 2. There, it was found that 
making inclusive cross sections exclusive by constructing explicit phase space subtractions 
(through the phase space mapping of PYTHIA8) will produce slightly harder partons in the 
core process (see Figure 11). This tendency persists after showering, as shown in Figure 
m for the UNLOPS case. Clearly, this is a non-negligible effect, although the differences 
are contained in the NLO scale variation band. We believe that using exclusive input is 
conceptually superior. However, this section is intended to give uncertainty estimates for 
NLO merged parton showers, and in particular to sketch merging scale uncertainty, and 
there is no reason to assume that the B- and the B-prescription differ in this respect. 
Using inclusive input, however, makes merging scale variations much simpler and quicker 
and avoids having to tamper with the internals of the POWHEG-BOX. Because of this speed 
factor, we chose to use inclusive input for the results of this publication. 

In the following, we will often include merging scale variations in the ratio plots. So 
that the plots become less cluttered, we will give the envelope of curves for merging scales 
between Pms = 15 GeV and Pms = 45 GeV as uncertainty band, rather than show the 
actual curves. 

Figure ^ shows that the transverse momentum of the hardest jet is heavily affected 
by NLO merging. This is due to the W+jet NLO calculation, as already seen in Figure |l]. 
The merging scale variations, as well as the //ij/^i?- variation for NLO merged results lie 
within the scale variation band of the NLO calculation, but the combined variation is not 
significantly smaller. The NLO merged predictions touch the upper limit of the NLO scale 
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Figure 5: Transverse momentum of the second hardest jet, for W-boson production in pp collisions 
at EcM — 7000 GeV, when merging up to two additional partons at LO, and zero and one additional 
parton at NLO. Jets were defined with the /c^-algorithm, with k±,rnin = 10 GeV. Multi-parton 
interactions and hadronisation were excluded. The lower inset shows the deviation from PYTHIA8. 
Left panel: Results of NL^. Right panel: Results of UNLOPS. 

variation band, because of the use of inclusive NLO input, as discussed earlier. Merging 
scale variations alone are minor. 

Merging scale uncertainties are also small in Figure |5|, which shows the transverse 
momentum of the second hardest jet. It is particularly reassuring that even combined with 
/Xij//ii?-variation, the bands are smaller than in CKKW-L and UMEPS [6] 

We would like to conclude this section by noting that differences between UNLOPS and 
NL^ are hardly noticeable for the displayed observables. This is true for all observables we 
have investigated in W-boson production, which can be interpreted to mean that the loga- 
rithmic improvements in UNLOPS do not result in major changes in W-boson production 
for the merging scale we have chosen. We anticipate larger effects once scale hierarchies 
become larger, i.e. if the merging scale is significantly decreased. For now, we will instead 
investigate if the introduction of a slightly larger mass scale and incoming gluons in the 
lowest order process leads to visible effects. 

4.2 H-boson production in gluon fusion 

This section is intended to demonstrate that the PYTHIA8 implementation is not specific 
to W-|-jets, and that different processes can be used to guide algorithmic choices. We 
have chosen to investigate Higgs-boson production in gluon fusion, mainly because of the 
presence of incoming gluons in the lowest order process and the very large NLO corrections. 

In Figure we compare the variation of NLO merged results with the scale variation 
in the H-|-jet NLO calculation of the default POWHEG-BOX program. The transverse 
momentum spectrum of the hardest jet is softer in NLO results than in PYTHIA8. We 
found the same behaviour in tree-level merging as well. Interestingly, a similar effect was 
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Figure 6: Transverse momentum of the hardest jet, for H-boson production in pp coUisions at 
EcM = 7000 GeV, when merging up to two additional partons at LO, and zero and one additional 
parton at NLO. Jets were defined with the fc^-algorithm, with k±^min ~ 10 GcV. Multi-parton 
interactions and hadronisation were excluded. The lower inset shows the deviation from PYTHIA8. 
The band labelled "POWHEG W+jet" is given by the envelope of varying the renormalisation 
scale in the POWHEG-BOX program between ^Mz, . . . , 2Mz, and the factorisation scale between 

^Mh, . . . , 2Mh. The upper panels show the results for using the zero-jet if-factor (i.e. K = J^) 
throughout the NLO merging procedures. The lower panel show the result when not using any 
if-factor (i.e. K = 1). Left columns: Results of NL^. Right columns: Results of UNLOPS. 

observed in pure QCD dijet production, which might indicate that PYTHIA8 overestimates 
the hardness of radiation from initial state gluons. 

We consider Figure ^ a cautionary tale. Let us examine the the upper row first. The 
merging scale variation in p±i is very small. However, when including renormalisation- 
and factorisation-scale dependence, the uncertainty of the NLO merged results is larger 
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Figure 7: Comparison of inclusive cross sections for H-production in gluon fusion pp collisions 
at Ecu = 7000 GcV, for NL^- and UNLOPS merging, when merging up to two additional partons 
at LO, and zero and one additional parton at NLO. The results of the merging procedure (labelled 
Cmcrgod) are normalised to the zero-jet inclusive NLO cross section CTinciusivo- The coloured bands 
indicate the variation from choosing /x_f in [ (^Mh)^, . . . (2Mh)^ ] and ^j.r in [ (^Mz)^, . . . (2Mz)^ ] . 
The error bars represent only the statistical error on the merged cross section. The results of the 
NLO merging procedure are presented for two different factor treatments. The result of using no 
K-ia.ctoi is labelled "no K-factor". Usage of the zero-jet iiT- factor is indicated by the label "0-jet 
K-factor". The same variations of NL'^ and UNLOPS are plotted, although the variance in the 
UNLOPS results is hardly visible. 

than the variation in the H-|-jet NLO calculation. This is explained by our choice of K- 
factor. As discussed in appendix this rescaling affects only the "higher orders" , since 
the effect on n-jet observables is removed to O (a"^^ (/^ij))- Different choices lead to no 
visible effects in W-boson production, since e.g. a change of Kq = = 1.16 w 1 will only 
result in changing tree-level samples slightly, and p±i is dominated by the one-jet NLO 
contribution. However, this does not apply to H production in gluon fusion, where Kq > 2 
leads to a significantly larger two-jet tree-level contribution. Enhancing the two-jet tree- 
level contribution will make the leading-order scale variation of this sample more visible, 
thus leading to an overall larger variation^^. 

Imposing a leading-order scale uncertainty on NLO observables is very conservative, 
and it seems prudent avoid artificial increases due to JT-factors that rescale higher orders. 
The lower row of ^ shows the result of not using any JC-factors at all. The agreement 

^^The same might naively be true for the POWHEG result, since two-jet contributions in H-|-jet in 
POWHEG also carry a (much more complex, phase-space dependent) T^-factor K — This "one-jet K- 
factor" increases with increasing fip.R and counteracts the decrease in H-|- 2 jet cross section with increasing 
scales. This leads - among other improvements - to a small scale variation in the POWHEG-BOX calculation 
for H-l-jet. 
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Figure 8: Jet multiplicity and transverse momentum of the hardest jet in W-boson production, 
as measured by ATLAS [46]. MC results taken from the POWHEG-BOX program, with three dif- 
ferent renormalisation/factorisation scales. Effects of partons showers, multiple scatterings and 
hadronisation are included. 



with POWHEG-BOX is reassuring, and the scale variation is small. The NL'^ result however 
exhibits major merging scale variations, which are mainly induced by an increased cross 
section in the results for = 15 GeV. This unitarity violation was previously "masked" 
by a large i^- factor. 

To further illustrate the effect of including iC-factors, we show the merged prediction of 
the total inclusive Higgs cross section as a function of the merging scale in Figure |7[ In this 
figure, we divide the NLO merged results for scales hr/ fip by the input NLO cross section 
with the same nn/fiF choices. The ideal result should be unity, without scale uncertainties. 
This is true to a high degree for UNLOPS, which shows that the unitary nature of that 
method really works as expected. For NL^, however, we see that when using JC-factors we 
get a large scale variation with a non-negligible merging scale dependence. Removing the 
-R'-factors decreases the scale variations, but on the other hand increases significantly the 
merging scale dependence. 

The i^-factor dependence is a major uncertainty in the NLO merged results for Higgs 
production in gluon fusion. We would like to stress that the current publication is intended 
as a technical summary, and not aimed at making binding predictions. Rather, we will use 
this as guidance when improving the implementation further. 

4.3 W-boson production compared to data 

In this section, we would like to show NLO merged predictions in comparison to data. 
We would like to point out that we have fixed as(Mz) in the PS to Qs(Mz) = 0.118, and 
use CTEQ6M parton distributions throughout. Please consult appendix ^ for a discussion 
of multiparton interactions. This means that the results do not correspond to a tuned 
version of the PYTHIA8 shower. Conclusive results can of course only be presented after 
the uncertainty of PS tuning has been assessed. 
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Figure 9: Jet multiplicity in W-boson production, as measured by ATLAS [46]. The MC results 
were obtained by merging up to two additional partons at LO, and zero and one additional par- 
ton at NLO. MC results are shown for three different merging scales (top panels) and for three 
different renormalisation/factorisation scales (bottom panels). Effects of multiple scatterings and 
hadronisation are included. Left panels: Results of NL"^. Right panels: Results of UNLOPS. 



In figure we show that the jet multiplicity is well under control in NLO merged 
predictions. The left panel of Figure ^ shows that, as expected, it is not possible to 
describe the number of zero-jet events with a W+jet NLO calculation. This is of course 
exactly the strength of merged calculations: Observables with different jet multiplicities 
can be described in a single inclusive sample. 

The transverse momentum of the hardest jet in association with a W-boson is shown 
in figure |l^ and the right panel of Figure ^. It is clear that the NLO merged results do 
not agree with data. We have chosen this particular observable because it our exhibits 
the most unsatisfactory description of data that we have encountered while testing our 
NLO merging methods. The reason for this disagreement is multifold. First, we have 
already mentioned that correcting for inclusive NLO input produces harder p±i tails. The 
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Figure 10: Transverse momentum of the hardest jet in W-boson production, as measured by 
ATLAS [46]. The MC results were obtained by merging up to two additional partons at LO, and 
zero and one additional parton at NLO. MC results are shown for three different merging scales 
(top panels) and for three different renormalisation/factorisation scales (bottom panels). Effects 
of multiple scatterings and hadronisation are included. Left panels: Results of NL^. Right panels: 
Results of UNLOPS. 



two-jet sample will eventually dominate the tail. We have chosen to rescale the two-jet 
contribution with a X-factor above unity. It could also be argued that the POWHEG-BOX 
result in Figure ^ has slight tendency to overshoot. This might indicate that some part 
of the "giant i^-factor effect" due to enhancements of O In ^|t"^ is developing in the 
W-|-jet NLO calculation of p_Li because of soft /collinear W-bosons. The last two points 
are correlated, since two-jet configurations have a major impact on the p^i-dependence of 
the NLO result, and increasing the two-jet contribution can enhance the visibility of giant 
i^T-f actors. 

The NL^ and UNLOPS descriptions of data exhibit high similarity. We have already 
noted the semblance of both methods in section 4.1 . This observation is specific to W-boson 



production, and does not hold for other processes, as for instance illustrated in section 4.2 
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5. Discussion and conclusions 



In this article, we have presented two new methods for combining multiple next-to-leading 
order calculations consistently with the PYTHIA8 parton shower. The NL^ method is a 
generalisation of the CKKW-L scheme, while the UNLOPS prescription accomplishes the 
same for UMEPS. Both methods achieve a description of zero-, one-, . . . , n-jet observables 
simultaneously at NLO accuracy in one inclusive sample, provided input event files at NLO 
accuracy for up to n additional jets are supplied. We would like to point the interested 
reader to appendix D.l, in which we argue out that it is feasible to extend the UNLOPS 
method to a NNLO matching scheme. 

Two distinct NLO merging schemes were presented to estimate the magnitude of issues 
related to sub-leading logarithmic enhancements. Although the UNLOPS method can be 
considered theoretically preferable, no large differences between NL^ and UNLOPS have 
been observed when merging multiple NLO calculations for W-boson production in associ- 
ation with jets. This leads us to conclude that for the observables that were investigated, 
and the merging scale values that were used, sub-leading logarithmic enhancements are 
sub-dominant. For H— boson production, differences are visible, with UNLOPS delivering 
a more reliable solution. 

This article is intended to give a comprehensive description of the choices that can 
be made in deriving and implementing an NLO merging method. We hope that this 
publication provides enough information about the actual implementation to allow the 
reader to form clear judgements of the rather intricate details. We have tried to remain as 
general as possible in our choice of inputs. It has been shown that different inputs can, due 
to mismatches in phase space mappings, have visible, systematic effects. When confronted 
with such effects, it is clearly preferable to reach an agreement over inputs, and we hope 
that the current publication can contribute to a discussion. 

We have shown that the merging scale dependence in W-l-jets is small, and contained 
in the scale variation band of the W+jet NLO calculation. This also means that the 
description of data is governed by the input NLO calculation. 

The merging scale dependence in Higgs-boson production in gluon fusion is very small. 
In this case, we highlighted that the dominant uncertainty of the algorithms is given by 
the choice of the if-factor rescaling higher orders - which is beyond the control of the NLO 
calculation. This is a manifestation of the magnitude of i^-factors in gluon fusion and the 
scale variation of the cross sections. We would like to stress that this uncertainty is present 
because we try to be as general as possible, and that the introduction of K-factors does 
in principle not jeopardise NLO accuracy, or degrade the PS approximation. However, if 
-ff-factors are not necessary and instead produce large variations, the removal of X-factors 
should be considered. 

Although we have presented some comparisons to data in this article, we do not attempt 
to make any definite predictions. To do this, a further investigation of the uncertainties 
has to be performed - a task we will return to in future publications. We end this article 
by listing the main issues that need to be addressed. 

Our methods require events generated according to the exclusive NLO cross section. 
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There are currently no standard programs that wih produce such events, and instead we 
have used inclusive NLO cross sections and subtracted explicit counter events by integrating 
tree-level matrix element events over the radiative phase space, using the mapping of the 
Pythias parton shower. We have also "hacked" POWHEG-BOX to directly produce the 
exclusive cross section event, and have found some differences, due to the different phase 
space mapping used there. Modifying the internals of other programs is, of course, not a 
viable long-term solution, and we hope that the introduction of our algorithm may inspire 
authors of NLO matrix element generators to include the generation of exclusive cross 
sections as an option in their programs. 

We have allowed the use of iC-factors in the underlying tree-level merging in the hope 
that the inclusion of NLO corrections will then lead to less merging scale variations. Al- 
though this can be done without modifying the formal accuracy of our methods, we see 
clear differences compared to the case where -fC-actors are omitted, in the case of Higgs- 
production, where these i^-factors are large. We find indications of reduced factorisation 
and renormalisation scale uncertainties in the absence of iiT- factor, but also note larger 
merging scale variations in the NL^ case. This needs to be investigated further. Other 
options, e.g. including multiplicity-dependent iT-factors, should also be considered to un- 
derstand uncertainties. 



At very high transverse momenta we expect logarithms of the form In j to arise, 
resulting in so-called "giant i^-factors" [36, 37]. These logarithms can in principle be 
resummed to all orders, and an inclusion of such resummation is planned for the parton 
shower in PYTHIAS. We are confident that our methods can be extended also to deal with 
this full electro- weak shower, but meanwhile we need to understand better the uncertainties 
arising from these logarithms. 

Finally, before we can be confident enough to make precise predictions with our new 
methods, a re-tuning of the shower (including MPI) of PYTHIAS must be carried out. 
The currently available tunes have all been obtained without higher order matrix elements 
merging, and it is clear that some of the resulting parameters have been obtained from 
trying to fit distributions where we do not expect an uncorrected parton shower to do 
a reasonable job. In particular, this applies to the tuning of the scale factor in Os (see 
eq. ( |2.8| )) in the shower, and we expect this to change significantly when tuning the ME 
corrected shower. This will then also directly influence the MPI, which also need to be 
re-tuned. Needless to say, such a tuning as a major undertaking. 

Note added 

While finishing this article, it came to our attention that an approach that is similar to 
UNLOPS has been developed in parallel by Platzer [32]. Also, on the day of submission, 
we noted that Aioli et al. [47] presented their work on NLO-matching. 
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Appendices 

A. NLO prerequisites 

This appendix is intended to introduce the merging scale definition used throughout this 



article (see section A.l ), discuss the prerequisites on NLO input (section [A.2| ), introduce 
the notation we employ (section |A.3| ) and finally illustrate how the POWHEG-BOX program 
can be used to generate the input necessary for NLO merging (section |A.4| ). 

A.l The Pythia jet algorithm 

Throughout this paper, we use cuts on the minimal PYTHIA8 evolution p±^evoh to disentan- 
gle regions of phase space. Since p±^evoi defines a relative p± distance [48] , we think of p±^evoi 
as an inter-parton separation criterion. To avoid confusion with other p^— definitions, we 
will use the symbol p for p±^evoi- 

The phase space regions in which we believe fixed-order calculations to dominate is 
separated from the resummation region by a cut value Pms, defined in a parton separation 
t ~ minj/o}. This minimal separation is constructed by finding the minimal p for any triplet 
of partons e, r, s, where e is a final state "emitted" parton, r is a radiating parton and s is 
a spectator. All triplets, irrespectively of fiavour (or colour) constraints are included. In a 
dipole picture, the radiator r can be thought of as the dipole end whose momentum changes 
most when splitting the dipole (r', s') into two dipoles (r, e), (e, s) while s is the dipole end 
that absorbs the (small) recoil. The functional definition of this parton separation criterion 
is 

t = min [/0{ij,A;}] where i G {final state partons}, (A.l) 
and j G {final and initial state partons}, 
and k € {final and initial state partons} 

where the separation of i for a fixed triplet k) of partons with momenta Pi,Pj,Pk is 

Zijki^ — Zijk)Qij if the radiator j is a final state parton, and 



_ 2 
Pijk — P±,evol,ijk 



12 _ /'^ I ^ \2 ^ _ ^i,jk 



(1 — Zijk)Qfj if the radiator j is an initial state parton, and 

^ij — KPi Pj) 1 Zi]k — (p-+ps^)2 

The cut value Pms is called merging scale. If all pijk for a particular final state parton 
i are larger than Pms, we call this parton a resolved jet. Conversely, if any pijk is below 
Pms) we call i an unresolved jet. We say that a phase space point is in the matrix element 
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region if t > p^s, i-e- all minimal parton separations are larger than the cut. In other 
words, a phase space point is in the matrix element region if it only contains resolved jets. 
The parton shower region is disjoint: If any jet separation falls below Pms, be believe that 
parton shower resummation is appropriate. 



Using eq. ( A.l ) and eq. ( A. 2 ) as merging scale definition does not exactly correspond 
to separating the matrix element- and parton shower regions in p±^evoi- In parton shower 
algorithms, the resolution scale attributed to a state is given by the scale of the last 
splitting. This is just one number, since a splitting is generated by a winner-takes-it all 
strategy: If a splitting is chosen, all scales attributed to splittings of other partons are 
considered higher. Such a merging scale definition can only be constructed if we know (all) 
parton shower histories of an input event. 



For now, we use eq. ( A.l| ) and eq. ( A.2[ ) as merging scale definitions, and are content 



with the fact that Pms does still correspond to a single p_L^evo«— value. This means that 
vetoing shower emissions that would result in an additional resolved jet will not introduce 
no-emission probabilities above Pms- 

We have to point out that fixing the merging scale definition is necessary in the 
NLO merging methods illustrated in this article. Otherwise, it would be mandatory to 
reweight NLO corrections with no-emission probabilities for merging-scale-unordered emis- 
sions, which would fundamentally degrade the higher-order description we aim to achieve^^. 
One benefit of CKKW-L tree-level merging is that the method allows for a wide class of 
merging scale definitions. Because of the treatment of emissions that are unordered in 
the merging scale, however, the merging scale effectively has to define a hardness-measure, 
since otherwise, only small portions of phase space will be endowed with ME corrections [5]. 
Different choices of hardness definitions for different processes in CKKW-L can be helpful 
for efficiently correcting phase space. The current implementation of CKKW-L in PYTHIA8 
allows for both minj/o} and min{A:^} as merging scales. No major efficacy differences be- 
tween these merging scale definitions has been found so far, leading us to conclude that in 
practise, defining the merging scale in min{p} is reasonable. 

A. 2 Exclusive cross sections 

In this section, we would like to introduce the concepts of exclusive and inclusive NLO 
cross sections, and comment on how inclusive NLO cross sections can be made exclusive 
by the inclusion of a phase space subtraction sample. 

We think of matrix element merging as a two-step measurement. We first measure 
the number of resolved jets in the input event, by applying a cut. Then, we calculate an 
interesting observable on events that have been classified as n— jet events. To make the 
second step independent of the choices in the first measurement, we need to sum over all 
possible jet multiplicities. 

Throughout this publication, we will define jets by the PYTHIA8 evolution jet 



separation criterion, as discussed in appendix A.l. Resolved jets are defined as partons 



^See appendix |^ for details 
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Figure 11: Comparisons of two way of generating exclusive NLO cross sections. Bi events are 
calculated without changing the POWHEG-BOX program, while for Bi, we have explicitly introduced 
the necessary cuts in PoWHEG-Box. Left panel: Comparison for zero-jet exclusive NLO cross 
section, as function of the W-boson rapidity. Right panel: Comparison for one-jet exclusive NLO 
cross section, as function the kinematical transverse momentum of the parton. 



whose separation to any other partons in the state is above the value Pms- We call Pms the 
merging scale. 

The output of a tree-level calculations for k final partons can contain I = 0, . . . ,k 
partons with soft or collinear momenta. The result of the calculation will diverge as soon 
as any parton in the calculation becomes soft or collinear. We can remove these regions 
of phase space, if we enforce that the output only contains exactly k partons with jet 
separation above Pms, i-e. k resolved jets. If the jet definition is infrared and collinear safe, 
this jet cut will render tree-level calculations finite. 

The result of using a next-to-leading order calculation for k partons to describe an 
observable O can schematically be written as 

{O) = J 0{cl)k) daBorn + J d'^cj^ 0{^k) davirtual + J d'^+V 0{^k+l) dcjReal , (A.3) 

where / d'^cf) indicates an integration over the fc— parton phase space, daBorn is the tree- 
level cross section, davirtuai the virtual correction term and da Real the real emission part. 



Equation A.. 3 allows contributions of any number of / = 0, . . . , k+1 resolved (or unresolved) 
jets. Since the tree-level part still diverges if any of the k partons approach the soft and 
collinear regions, we require that 0{(j)k)daBorn always contains exactly k resolved jets, 
which immediately means having the same requirement in J^^O{(j)k)dayirtuai- This then 
means that J^^-,^ 0{(t)k+i)daReai has to be constrained, since otherwise, the NLO calculation 
could include real-emission corrections to non-existent "underlying" Born configurations 
with less than k resolved jets. The POWHEG method [17, 18] eliminates this issue by 



evaluating using one phase space point cp^ for tree-level and virtual parts, and (pn+i phase 
space points that can be projected exactly onto (p^^ in real-emission terms. We will assume 
that in neither of the terms in |A.3 , less than k resolved jets are included. Further, the 



observable O receives - through real corrections - contributions from A; -|- 1 (resolved or 
unresolved) jets. Measurements explicitly depending on the kinematics of A; -|- 1 resolved 
jets are only accurate to tree-level approximation, while contributions for an unresolved 
additional jet are necessary to cancel divergences^^ 

In multi-jet merging, (k + 1)— jet contributions enter through explicitly adding a 
reweighted {k + 1)— jst sample. To merge multiple NLO calculations, we also need a clean 
cut, that makes the classification of the input in terms of jet multiplicities possible. We 
define that the "NLO part" of an A; —jet NLO calculation should contain k resolved jets, 
and at most one unresolved jet, while contributions with A: -|- 1 resolved jets are regarded 
leading-order parts 

(O) = j d^cf) 0{(l)k) I da Born + davirtual + j d<t) daReal^ (A.4) 

^ V ' 

NLO part 

PMS 



LO part 

We call the NLO part of such a calculation the exclusive NLO cross section. 

Merging schemes naturally act on exclusive cross sections. Technically, we assume 
that all unresolved real emission parts are projected onto n-parton phase space points. 
In POWHEG, this is facilitated by performing the d(j) integration in the real-emission term 
explicitly. If an NLO calculation would yield only phase space points with exactly n partons 
in n resolved jets, and weight these points with^^ 

^ rpMS 

Bn = daBorn + dcJvirtual + dcj) da Real , (A. 5) 

then this calculation could be used immediately for NLO merging. However, B„ might 
not be accessible without changing the NLO matrix element generator, since it might not 
possible to split the calculation into "NLO parts" and "LO parts" as desired. We instead 
choose to write 

B„ = daBorn + davirtual + I dcj) daReal - / # daReal (A. 6) 

J J pMS 

and to use the two samples 

Bn = daBorn + dayirtual + J dcp daReal (A. 7) 

B„+l^n = - dcf) daReal (A. 8) 



^^In numerical implementations, the singularities of virtual corrections and real-emission contributions 
are cancelled separately by regularisation terms. 

^*We give more precise definitions of the notation in the next section. 
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separately. We call B„ the inclusive NLO cross section, and refer to j\, B^+i^n as phase 
space subtraction term. Adding the inclusive NLO cross section and the subtraction term, 
we retrieve the exclusive NLO cross section. It is possible to formulate NLO merging acting 
on exclusive cross sections, or acting on inclusive cross sections and additional phase space 
subtraction samples. 

We would like to comment on our framework to generate exclusive cross sections from 
inclusive NLO input. It has already been demonstrated in [6] that extracting integrated 
states that had been constructed as intermediate states in the parton shower history does 
indeed give the expected results (see Figure 2 of [6]). However, these comparisons were per- 
formed between the shower approximation and the integrated matrix element. To correct 
for the use of an inclusive NLO cross section, we have to use integrated tree-level events 



as phase space subtraction for POWHEG-BOX events. Figure 11 shows how the exclusive 
cross section generated through explicitly changing the POWHEG-BOX program code com- 
pares to the exclusive cross section produced by a-posteriori phase space subtraction. The 
differences between the two prescriptions stems from a different phase space mapping in 
POWHEG-BOX and PYTHIA8. Such effects are beyond the accuracy of the NLO merging 
methods presented in this article. It should be noted that an a-posteriori phase space 
subtraction using the PYTHIA8 phase space mapping produces a harder p^-spectrum and 
less forward W-bosons. Ideally, we would like generate the exclusive NLO cross section 
with an NLO generator which allows for the necessary cuts, so that explicit phase space 
subtraction become unnecessary. For this publication however, we will use subtractions, 
since this allows us to perform merging scale variations without continuously having to 
re-generate LHEF output with the POWHEG-BOX program. 

A. 3 Notation 

Formulating NLO merging is unfortunately a fairly notation-heavy task. In this section, 
we would like to carefully introduce the symbols we will use throughout this article in a 
tabular style, for easy reference. Let us define the lingo 

(pn- Phase space point with n additional resolved jets. This 

means that can contain p lowest-multiplicity particles 
(e.g. e"*" and e~ for Drell-Yan production), and n additional 
partons. Each phase space point has a fixed momentum, 
flavour and colour configuration. 

Underlying configuration: Phase space point (pn which can be to constructed from 

(pn+i by removing one emission, meaning integrating over 
the one-particle phase space of the emission, and recombin- 
ing flavours and colours. We will use the terms underlying 
momentum configuration, underlying flavour conflguration 
and underlying colour conflguration when explicitly empha- 
sising one aspect of the underlying conflguration. 
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The radiative phase space, meaning that in 



(l>n+l ~ 4>n^rad 



(A.9) 



Resolved jet: 



$rad plays the role of the one-particle phase space of the 
additional parton, while (pn is the underlying configuration. 
Parton, for which all possible jet separations to other par- 
tons in this phase space point are larger than the value 



Unresolved jet 



Pms- 

Parton, for which at least one jet separations to one other 
partons in this phase space point is lower than the value 



Pms- 



We always assume that matrix elements for n outgoing partons arc only integrated over 
phase space regions with exactly n resolved jets, unless explicitly stated otherwise. In 
cases where we integrate over emissions, we assume that in the ME input, all partons were 
resolved jets, and that after the integration, all partons are resolved jets. 

Please note that the methods presented in this publication do not rely on a particular reg- 
ularisation scheme in the NLO calculation, as long as the dependence on the regularisation 
scheme is cancelled locally, i.e. in the weight of each phase space points separately. We 
will choose a rather symbolic notation for parts of NLO calculations, and hope that this 
will make the formulae in this article more accessible. In the following, we would like to 
introduce the shorthands: 

B„: Tree-level matrix element with n partons, i.e. 



Bn = f^ix+,IJ,F)fn {Xn,IJ'F)\Mn,o{lJ'F,IJ'R)\'^ , (A.IO) 



where the first subscript on |7Wn,o indicates the number of jets, 

while the second counts loop integrations. Readers more familiar with the 
notation of [18] should note 



For brevity, we always suppress flavour indication on B; 
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B, 



n+l\n- 



Js n- 



Sum of tree-level configurations with n + l partons, for which the underlying 
Born configuration is (pn- Loosely, we think of ^n+i\n the sum of matrix 
elements B„, multiplied by splitting kernels (with B„ being evaluated us- 
ing the underlying momentum, flavour and colour configuration, while the 
splitting kernels depend on the radiative phase space). Translating to the 
notation of [18], this means 



n+l\n 



E 

ar^{ar\fb} 



n+l) 



We choose the more symbolic indexing with n + V\n, rather than the more 
rigorous ar-notation of [18], in order to not obfuscate formulae with details 
that are not essential for the discussion. 

Sum of tree-level cross sections with n resolved jets in the input ME events, 
after integration over the phase space oin — m partons. Symbolically: 



B. 



E 

n 



B 



n\n—m 



(A.ll) 



We think of being produced from ^n-m by m consecutive splittings. The 
index s denotes that the integration is accomplished by explicitly removing 
m partons from the n— parton phase space, meaning that we substitute the 
state S'+n by the state 5*+ 

n— m? introduced in [6]. The symbol J"^ B^_>m 
also indicates that the state S+m only has resolved jets, and that possibly 
more than one integrations had to be performed in case all of the states 
S+n-i, • • • 5 S+m+i contained unresolved jets. 

Sum of tree-level cross sections with n resolved jets in the input ME events, 
after integration over the phase space of n—m partons. However, in contrast 
to the symbol B^^m, we explicitly require the states S'+n-i to contain n — 
1 resolved jets, and still perform this second integration. This is indicated 
by the upward-pointing arrow. All further integrations only have to be 
performed because the states S+n-2, ■ ■ ■ , S+m+i contained unresolved jets. 

Bg^Q for example means that we first replace the state S+s by 5+2 (with 
two resolved jets), then demand another integration, giving S'+i. Then, we 
find that S'_|_i contains an unresolved jet, so that we integrate once more. 
The last step would not be necessary for the term B^ 



Virtual correction matrix element with n partons above Pms- 



Vn = fniXn,IJ'F)fn , I^f) \Mn,l {iJ-F.IJ'R) 



(A.12) 



We assume that all ultraviolet divergences have already been removed. 
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Dn+i|n- Sum of infrared regularisation terms with n partons above Pms- As above, 
we indicate these terms can be projected onto underlying Born configura- 
tions by the index n + l\n. For simpUcity, we may think of these regulators 
as Catani-Seymour dipoles, and very symbolically put 

i'j'k 

(g) Dij^i/ji^. ($rad) ) 

where i and j are partons of the underlying configuration 0„, while i' , j' and 
k are partons of (pn+i- As long as all dependence on the regularisation is 
contained in the inclusive (exclusive) NLO cross sections, our method will 
not depend on the actual form of these terms - all numerical NLO subtrac- 
tion schemes are equally valid. The notation n-|- l|n is to be understood in 
the same way as for B„_|_x|„. 
In+i|n- Sum of integrated infrared regularisation terms with n partons above /9ms- 
Remainders due to initial state partons being collinear with identified initial 
hadrons are included in I^+ii^j. As above, we indicate that the terms in 
ln+i\n can be projected unto underlying configurations by the index n + l\n. 
Schematically 

n+l\n 

(A. 14) 

In this case, the one-parton phase space integration is commonly performed 
analytically. The integration here covers the complete radiative phase space. 
Integration variables can change for different types of dipoles Dij^i'j'^. 



- 42 - 



B„: Inclusive NLO weight of ra— parton phase space points. 

Bn = B„ + V„ + I„+l|n 

+ J d^rad (Bn+l|n - Dn+l|n) • (A.15) 

Note that we will assume that this gives a NLO weight for the phase space 
point 4>n, meaning that we have to evaluate ^n+i\n ^n+i\n with rather 
than (fn- This is the standard procedure in the POWHEG and MC@NLO 
methods. The integration over $rad covers the full radiative phase space. 
Real emission terms can give contributions with an additional resolved jet, 
which are here included into the weight of n-jet phase space points 
As discussed above, the description of an observable at NLO receives n-jet 
and an + 1-jct contributions, so that projecting all n+ 1-jet configurations 
onto (pn (and leaving no n + 1 events) seems problematic. However, in 
matrix-element merging, n + 1-jet events will be included though the next- 
higher multiplicity sample. We then need to ensure that the contribution of 
resolved n + 1-jet events to the cross section is not double-counted. This is 
solved by the introduction of exclusive NLO jet cross sections. The inclusive 
cross section is closely related to the definition of B in the POWHEG method. 
Indeed, if no Sudakov factors are applied in the weight of Born-type phase 
space points in POWHEG, and all radiative events are projected unto Born 
configurations, this exactly produces our definition of B„. 
B„: Exclusive NLO weight of n— parton phase space points, 

Bn = B„ + V„ -|- I„+i|n 

/Pms 
d^vad (B„+i|„ - D„+i|„) 

= B„ — / (i$].adB„+i|„ 

= B„-^B„+i_^„ . (A. 16) 

It is clear from the last equality that an exclusive NLO n— jet cross section 
can be constructed from the inclusive case by explicitly subtracting the 
phase space points with an additional resolved jet. 
Bn-^m- Inclusive NLO cross section with n resolved jets in the input ME events, 
after integration over the phase space of n—m partons. The symbol B„_^m 
as always also indicates that more than one integrations had to be performed 
because all of the states S+n-i, ■ ■ ■ , >S'+to+i contained unresolved jets. 
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Jin^jn- Exclusive NLO cross sections with n resolved jets in the input ME events, 
after integration over the phase space of n — TTi partons. The symbol B^^— 
as always also indicates that more than one integrations had to be performed 
because all of the states S+^-i, • • • , S+m+i contained unresolved jets. 

We will use several different event samples as input for multi-jet merging. An event consists 
of a phase space point with an associated weight, and can thus be considered completely 
differential. Only if necessary will we talk about predictions for an observable. We will use 
"cross section" and "event" interchangeably, and also make little distinction between the 
terms "(phase space) weight" and "matrix element". 

Exclusive cross sections are the basic building blocks needed for multi-jet merging. 
Tree-level merging uses phase space points weighted with exclusive tree-level matrix ele- 
ments B„ as input. The NLO multi-jet merging prescriptions advocated in this publication 
analogously require phase space points weighted by exclusive NLO weights as input. If no 
exclusive calculation is available, it is possible to extend the algorithm to include the ex- 
plicit subtraction of eq. ( A.16| ). To make formulae for NLO merging a bit more transparent, 



let us introduce the short-hands 

B„: UMEPS-reweighted tree-level cross sections with n resolved jets in the 

input ME events. 

Bn^rn- UMEPS-reweighted tree-level cross sections with n resolved jets in the 
input ME events, after integration over the phase space oin — m partons. 
The symbol B 

n^m also indicates that more than one integrations had 
to be performed because all of the states 5+^-1, • • • , 5'+m+i contained 
unresolved jets. 

[M-a b- Contribution A, with terms of powers a" and removed. The resulting 

terms are calculated with fixed scales hr and fip. 
[A]^ ^: Contribution A, with only terms of power and retained, calculated 

with fixed scales fiR and fip. 



The last two short-hands are particularly useful when trying to summarise terms in the 
expansion of the tree-level merging weights. For example, the sum of the second and third 
term in curly brackets in 



B2W2 

is given by eq. (B.46) below. 



B 



-2,3 



2 \^W2 



W2 



A. 4 Povi^heg-Box usage 

This section is intended to give guidelines on how to use the POWHEG-BOX program [19] 
in order to produce the inclusive NLO cross sections [41-44] needed for NLO merging in 
Pythias. Ideally, this should suffice as tutorial on how to set up the desired POWHEG-BOX 
outputs. We rely on knowledge on POWHEG-BOX input manipulations. With new versions 
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of POWHEG-BOX, the names of the inputs might change, so that the settings advocated 
here come without guarantees. 

For a given process with n partons in the underlying Born configuration, POWHEG-BOX 
by default generates output Les Houches events with n— and (n + 1)— parton kinematics, 
n— parton phase space points are weighted with the matrix element weights 

Bn A (p_L,mm) 

= \ B„ (<^„) + V„ + (A.17) 



+ 



J i^rad [B„+i|„ (0n (l^rad) - ^n+l\n {<t>n 0rad)] i A (p±,mm) , (A.18) 



where the integration / o?$rad contains the complete radiative phase space, and the Sudakov 
factor is given by 

A(PJ = <^-p[-lf,^ B„(0„) 1 ■ <^-''" 

{n + 1)— parton phase space points are weighted with 

B„ A ip^) ^"t!'",iV'^ e {p^ - . (A.20) 

Bn [<Pn) 

The program decides if a radiative (i.e. (n+ 1)— parton) phase space point is generated by 
comparing the p± of the proposed configuration (pn+i against p±^rnin- No radiative events 
are produced if p±,mm is set to the kinematical limit. Furthermore, we have 

{^(^^^.--)W„->oo = l- (A-21) 
Thus, using p±^rnin — >■ oo, we find 

{B„ A {P±,min)}p^^^.^^^ = Bn (</>n) + (<^„) + {<Pn) (A.22) 

+ J d^rad [B„+i|n {(pn ^rad) — D„+i|„ {(l)n (prad)] 

This is exactly the inclusive NLO cross section we need to perform NLO merging. The 
NLO merging prescription will include (n + 1)— parton configurations in a CKKW-L style. 

In the POWHEG-BOX program, the parameter p±^min can be set by changing the input 
variable ptsqmin. For example assigning 

ptsqmin = ldl5 (A.23) 

will ensure that for LHC energies, the output events of POWHEG-BOX will contain only 
n— parton kinematics, weighted with the desired inclusive NLO cross section. Setting 
ptsqmin will produce only n— parton kinematics only if every (n + 1) — parton phase space 
has an underlying Born configuration. This is for example not true for cc — )• udW~ 
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scattering via an s— channel gluon. Since such processes do not constitute corrections 
to any lower-order process, we regard these as true leading-order parts, and (internally) 
neglect these configuration in the POWHEG-BOX output. They will be added by including 
such configurations with incomplete parton shower histories through the treatment of tree- 
level matrix elements. Numerically, the treatment of incomplete states does not have any 
impact. 

The desired POWHEG-BOX output should be generated with fixed factorisation- and 
renormalisation scale. To be completely certain that this is the case, set 

runningscale (A. 24) 

runningscales (A. 25) 

btlscalereal 1 (A.26) 

btlscalect 1 (A.27) 

ckkwscalup (A.28) 



The POWHEG-BOX program would, upon making the assignment A.23| , attribute a SCALUP 



value of SCALUP = ^ptsqmin to the output LH events. This number would normally be 
read by PYTHIA8 and used as factorisation scale in the construction of overestimates for 
initial state splittings. For our purposes, the true value of fip will be an input for PYTHIA8, 
so that the correct choice can be used internally. 

After these settings, the POWHEG-BOX output file can be used for NLO merging in 
Pythias. We will include a detailed documentation of this procedure, and a manual how 
the schemes presented in this publication can be used, in the online documentation of an 
upcoming PYTHIAS release. 

B. Generation of weights 

The aim of this appendix is to provide a complete description of the O (og (^_r)) —terms 
needed to implement the NL^ and UNLOPS schemes. This task is split into subsections 
containing the expansion of factors appearing in the weight Wn 



2 = 1 



^S+niXn,Pn,PMs) , ■ (B.l) 



The weight applied in UMEPS differs in that the last no-emission probability 
x-m Pm Pms) is not included. Note that we have kept the parameter bi discussed 
in The generation of the tree- level weights for CKKW-L is discussed in [5], and the 
UMEPS case is treated in [6]. The iiT-factor is generated by dividing the (integrated) 
inclusive NLO zero-jet cross section by the leading-order result 



-46- 



It is in principle possible to rescale the tree-level weights for n partons by X-factors depend- 
ing on the jet multiplicity. Multiplicity-dependent X-factors will lead to different rescaling 
of PS higher orders, since the O {ij-r)) —contribution from multiplying i^-factors will 
be removed. This means that different ii'- factor choices give changes beyond the accuracy 
of the methods. It is interesting to observe that MC@NLO and POWHEG also differ in the 
X-factor applied to real-emission events: While the radiative events are not rescaled in 
MC@NLO, POWHEG includes a phase-space dependent iC-factor (see eq. ( |A.19[ )). Though 
formally sub-leading, this difference can be large, particularly if the NLO result is sig- 
nificantly higher than the Born approximation, e.g. in Higgs-boson production in gluon 
fusion. An argument against n— dependent JC-factors is that rescaling every jet multiplic- 
ity in CKKW-L by different numbers will result in increased merging scale dependencies. 
In Pythias, we include the possibility for having n— dependent X-factors 

Kn = . (B.3) 

■J PMS 

These will then be calculated by dividing the sums of inclusive NLO and LO cross section 
weights (with no extra reweighting) of events above the Pms cut. 

Below, we will give a detailed expansion of all factors in the tree-level weights, which de- 
pend on as, directly and indirectly through the PDF-ratios. Since we have demanded that 
the input cross sections be calculated with fixed fip and fiR, we will only keep contributions 
of powers al, and fixed scales. Incoming particles with positive (negative) momentum com- 
ponent pz will be indicated by a superscript + (— ). Final state partons will be enumerated 
by the superscript k. We will make use of the short-hands 

j&{q,q,g} 

ie{(?,9,g} 

where Pji are the unregularised Altarelli-Parisi splitting kernel for an initial state parton 
changing from i to j (by backward evolution) , and Pij are the plus-prescription-regularised 
counterparts for forward evolution. For final state splittings, we will write 

If the PS step from S'+j to Sj^i was final state 
Pji{y) radiation off leg k, with a final state recoiler. 

P^iiv^p) = * BS step from 5*+^ to 5*+^ was final state 

Pji{y) min i 1, \ ^ — r- \ radiation off leg k, and involved the incoming 

parton zb as recoiler. 

(B.7) 
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The factor mm{l, ^f^{^, p)/x^ f^{x^ , p)} is introduced on purely technical grounds, 
because the overestimate of final state radiation in PYTHIA8 does not include PDF fac- 
tors, and violations of the overestimate need to be avoided. We split the expansion of 
eq. (B.l) into subsections containing detailed expansions of each factor. At the end of each 



subsection, we will give a description of how the necessary terms are generated in PYTHIA8. 
B.l Expansion of iT-factors and ag— ratios 



The weight Wn in eq. ( [B.l[ ) contains the factors 

asikPi) 



Note that we have kept the parameters bi G {bijbp} stemming from different 
as(Mz)— values in parton shower and fixed-order calculation. The factors have simple 
Og— expansions 

K = 1 + a,{pR)ki + 0{al{pR)) (B.8) 
as(6i/0i) = «s(/ii?)|l + ^asMln(^^^| + 0{al{pR)) (B.9) 

where /3o = 11 — Multiplying these series, we get the expansion of the product of 

X-factors and Og— ratios 

n = 1 + aMki + E o^s{m)T- 1^ (t^) + 0("s (/^i?)) (B.IO) 

i=i '^s\^^R) ^'^ \OtPi/ 

We generate the /ci— term by using ki = K — 1. The sum is generated by stepping through 
the chosen PS history, and adding, for each nodal state S+j, the logarithmic terms, evalu- 
ated at the reconstructed splitting scale pi. This of course means that we have to construct 
and choose a parton shower history first. 

B.2 Expansion of ratios of parton distributions 

The expansion of the PDF ratios 

+ + + _ _ _ 

i=l 

is a bit involved. To derive an expansion, we will infer the DGLAP equation 

Jpi-i-Sp P Jxf y 
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where we have used the notation of eq. ( B.4 ). So far, no approximation to DGLAP scale 
dependence is made. This equation has contributions from all orders in as, and should 
be regarded as an expansion in the difference of scales 5p, rather than an expansion in 
as- Approximating fi{^^p) the sum of products of PDFs and leading-order or next- 
to-leading order splitting kernels (indicated by superscripts (0) and (1) respectively), we 
find 



(B.13) 



Shifting the scale in as to ^r, and the scale of parton distributions to fip, this gives 



+ 



+ 



+ 



+ 



«s(Mr) 



Pi-i 



p A± y ^ 



(B.14) 



P'-' dp dy 



dp' 



PF 



dy' 7±,(o). 



t/y y 



Rj 



2tt 
27r 



Pi-i 



PF 



p Jx^ y ^ 



+ 0{al{pR)) 



where f^'^^^ is a convolution of parton densities and leading-order DGLAP splitting kernels 
(see eq. (p.4D), and f^'^^^ a convolution of f^'^^^ and splitting kernels. 



The weight Wn contains ratios of parton distributions. Since some of these ratios are 



the result of rescaling to the lowest-order cross section (see discussion after 2.2), it might 



well be necessary to divide next-to- leading order PDFs. Luckily, eq. (B.14) ensures that if 
we are only interested in the expansion up to O (q^ (pr)), we can safely ignore difficulties 
relating to NLO splitting kernels. Using 



P^-^ dp 

PF P 



In 



Pi-i 

PF 



(B.15) 



and restricting ourselves to O (a\ (pr)), we arrive at 



f^{xf,p,.i) = ft{xf,pF) + 



as(/i 



R. 



27r 



In (^) [[ ^fi"\f,PF) + 0{al{pR)) (B.16) 



From now on, we will drop the superscript (0). With eq. ( [B.16| ), the expansion of a ratio 
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of parton distributions is given by 



xf_if^iixf_^, Pi) 27r 



R) 1 Pi-i 



dyxtJtii^y^:^ 

y xf_j^_^{xf_^, fxp) 



+ Oial{pR)) 



(B.17) 



With this, we can write the expansion of the product of PDF ratios to O (hr)) as 



ifn (Xn ,Pn) TT ^'i- 1 Ji- 1 l^i- 1 '/ 
fn (Xn ,IJ,f) J- J- X+_^f+_^{xl_^ 



1=1 



.Pi) a:j_i/j_i(3;j_i,pi) 



= 1 + 



27r 



In 



{^} 



dyX+f+{^,PF) 

y xt fn {xt , IJLf) 



XP'F) I y Xn fn {Xn , fip) 



+ 



E 



1=1 



In 



i\ I dy xt_Jt^{^,PF) 



+ 



E 



1=1 



In 



dyXi_Ji_^{^,HF) 
y x-_^f~_^{xr_^, hf) 



+ 0{al{iiR)) 



(B.18) 



These integrals can be calculated by explicit numerical integration. Remember that / has 
been defined with regularised splitting kernels [49] 



^qq(^) = (13^ + ^'^^^^^ = ^^^^'^ 

Pg^(^) = C^l + ilZfl! = Pqq(l-Z) = Pgq(z) 



(B.19) 
(B.20) 



Pgg(^) = 2Ca 



1 -z 



+ z{l - z) 



1 



Pqg(z) = Tr [z^ + (1 



+ -[llCA-^nfTR]5{l-z) (B.21) 

(B.22) 
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By using these functions explicitly, and inserting 

dy 

Jo 



1-2/ 



ln(l — X 



(B.23) 



we find that in the case that i — 1 is a quark or antiquark, the integral in eq. (|B.17D becomes 



2tt 



Pi 



dsifJ-R) Pi-i\ I dy ^i^ifi^i{^,PF) 



Ois{PR) I Pi-1 



27r 



dy 



+ / dy 



Tniy^ + (1 - yf) 



y My y ,Pi 



Xi-lf\{Xi^l,Pi) 



2Cf 



y y 'P' 



+ 2Ci.ln(l-x,_i) + -Cf 



(B.24) 



To arrive at this form, we have used that the sum over all possible flavours in reduces 
to two terms, since the antiquark (quark) can only be produced by the evolution of a gluon 
or an antiquark (quark). If the flavour index i — 1 indicates a gluon, the integral reads 

1 ^ 

dy xtJtiC-^,PF) 



asiPR) ) Pi-i 



2tt 



Pi 



y xf_Jt-iixf_^^ pf) 



asipn) ( Pi-i 



dy 

i-y 



2Ca y 



y ■'sl y 



Xi-lfg{Xi-l,Pi) 



2Ca 



+ / dy 2Ca 



i-y 



+ y(i-y) 



y ■'sy y 

Xi-lfg{Xi^l, Pi) 



Xi-l 

1 



+ dyCF 



i + (i-y)' 



Xi_l/g(Xi_l,Pi) Xi_l/g(Xi_l,/>i) 



y q \ y 



+ 2Ca ln(l - Xi^i) + - [IICa - AufTR ] 



(B.25) 



- 51 - 



The sum over all flavours is reduces to three terms - the evolution of a gluon, a quark 



and an antiquark. Since all remaining integrands in eqs. B.24 and B.25 involve PDFs, we 
need to perform the integrals numerically. This numerical integration will be performed 
by Monte-Carlo integration, although we have also implemented the Gaussian adaptive 
quadrature method for this task. We checked that compared to performing a many-point 
numerical integration, it is sufficient to only pick a single integration point for the Monte- 
Carlo integral evaluation, already when averaging over a small number of events (0(50)). 
Thus, we choose the less time-consuming Monte-Carlo method as default. For this, we roll 
a uniformly distributed random number r G [0, 1], and pick the integration variable as 



Vr, 



(B.26) 



if i — 1 is a gluon 
+ r • (1 — x) otherwise 

The full weight in eq. ( |B.18| ) is generated by stepping through the chosen PS history, and 
adding, for each reconstructed state 5+1-1 — 1 < n), the bounded integral for both 

incoming partons, multiplied by the logarithm of the reconstructed scales and pi. The 
values Xi-i, and pi are easily accessible since the parton shower history contains a 
complete sequence of fully reconstructed states S^q, . . . , For the ME state we 
add the x„— bounded integral, multiplied by the logarithm of pn and pp, for both incoming 
partons. 

B.3 Expansion of no-emission probabilities 

The factors Hs'^. j (xj_i, p^) in eq. (|B.1| ) symbolise the probability of a state 
to evolve from scale pi-i to pi without resolving emissions. This means that 
^s+i-ii^i-i, Pi-i, Pi) is a product of no-emission probabilities for each dipole in the state: 

n5+,_i(xi_i,Pi-i,Pi) = U+^^_^{xf_^, pi_i, Pi) U-^^_^{xr_^,pi_i,pi) (B.27) 

n ^s+.-^i4-vPi-uPi) 

fcsjfinal partons} 

Using the notation of eq. ( |B.5| ) and eqs. ( |B.£| ), we can write 

/■ Pi-l ^ 

dp I dy as{bip) xf-ifi^ii-Y', p) 
P J y xf_jt-ii^f-v P) 



n 



(B.28) 



and 



/- Pi-l 



n|+,_i(x^i,Pi-i,Pi) = exp< 
From eq. (B.29), we find 



dp / dy as{bpp) 



ply 2tt 

Pi Qp 



PU{y.p) 



(B.29) 



Pi~i 



ni,_,(xf_i,P.-i,P.) = l- / ^ / ^-^'^^PUy,p) + 0{al{p)) . (B.30) 
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All NLO contributions included in NLO merging schemes are generated with fixed scales 
and fip. For a NLO-accurate description, we need to remove the approximate shower 
contributions for exactly these scales. Otherwise, e.g. if we choose to remove the sec- 
ond term on the right-hand side of eq. (B.30), we will remove O (a^ (fJ-R)) terms as well, 
thus degrading the description of higher orders. To remove only precisely those parton 
shower terms that have corresponding contributions in the NLO calculation, we extract an 
O (ag (a^-R)) —expansion from eq. ( 6.301 ): 

Pi-l 



27r 



dp j dy 

P J y 



Pl^{y,PF) + 0{al{^,R)) ,(B.31) 



where we have used that the difference between OsI^fA*/?) and as{pR) is of O {pr})- 
The PDFs in appearing in P^_i{y), for final state radiation with an initial state recoiler, 
should be evaluated a,t pp. Expanding the exponential in eq. (|B.2^ ), we find 



pi-i 



dp J dy as{bip) xf_if^i{^^, p) 



P 



y 2vr xtjtii^tvP) 



+ 0{al{p)) , (B.32) 



which can be expanded further to give 

Pi-i 



n 



S+i_S^f-l^ Pi-l^ pi) = 1 



ols{pr) I dp I dyXi_if., 



± 7± ^^f-i 



• Pf) 



27r 



+ 0{al{pB)) .(B.33) 



P J y xf_Jt^^{xf_^, pf) 
Pi fii 

Here, we have used the result eq. ( B.l(j ), i.e. that if we are interested in the O (a^ (pr)) 
parts in eq. ( |B.33 ), we can safely evaluate the ratio of parton distributions at pp, rather 
than p. 

When expanding the CKKW-L weight Wn, we also need to discuss the "last" no- 
emission probability 

Us+„{Xn,pn,pMs) ='n^^Jx^,Pn,pMs) 'ns^Jx-,Pn,pMs) JJ ^S+Jxt^Pn,pMs) (B.34) 

fcG {final partons} 

with 



^s+S^t^Pn,PMs) = exp < 




dp I dy as{bip) ^nfni^^P) 



Pjy 27r Xnfn{Xn,p) 

, Pn+l 



e{t{S+n+i,p)- Pms) } (B.35) 



and 



r P' 



n|+„(x^,Pn,PMs) = exp < 



dp I dyasibpp) 




p)e p) - Pms) > (B.36) 
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The expansion of these terms carries an additional Q function 



(B.37) 



^^+„i^n ^ Pn, Pms) — 1 



Q^s(/^_r) I dp j dy ^nfni^j P-f) 

27r J Pj y Xnfn{Xn,fJ.F) 

Pn+l 



@ {t{S+n+l, P) - Pms) + ••• 



and 




^S+n(^nJ Pn, Pms) — 1 



P::{y,PF)@{t{S+n+i,p)-PMs) + ... (B.38) 



The definition of t is given in appendix A.l. Collecting all terms, we find that the expansion 
of the product of no-emission probabilities in the CKKW-L weight Wn is given by 



^S+A^n, Pn, Pms) Y[ (xi-i, Pi-i, pi) 



(B.39) 



2=1 



Pi-1 



Eh: 



a. 



{Pr) dp dy4-iftiC-^^PF) 



2tx 



± 7± 

y_ 

p J y xf_jtii^f-vPF) 

Pi 0,1 

pi-1 



OisiPR) dp dy 



2tt 



p y 



Pi-iiy.PF) 



Pi Qf 



pn 



sr^asipR) I dp I dyXnfni^,PF) 



2ir 



P I y xtfn{Xn,PF) 

Pn+l 
Pn 



E 



as{pR) dp dy~k 

J pJ y 

Pn+l Clp 



27r 



P^{y,PF)Q{t[S+n+i.p)- Pms) + 0{al{pR)) . 



Although this looks fairly complicated, it is in fact easily generated. It is useful to remember 
that if the probability for n incidents (e.g. nuclear decays) is given by 



1 



(B.40) 



then the average number of incidents is 

oo oo oo ^ oo 

=1 ^ n=l ^ '' 



n=0 n=l 
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The probability P„ has the same form as the probabihty for n emissions in a parton shower, 
and the O (og (/^i?)) —terms in eq. ( [B.39| ) can be identified with the exponents x. That 
means we can calculate the O (og (fJ-R)) —terms by generating, with fixed hr and /i^, an 
average number of emissions. In final state radiation off final parton k for example had 
been generated with fixed scales, we can write 



emissions between pi—i and pi 
Pi-i 




Ply 



exp 



pi-i 




p y 



p y 



(B.42) 



The average number of emissions is additive. This means that instead of averaging over 
emissions from each leg separately, we can directly average all emissions: We can simply 
start the shower off state 5'+i_i at Pi^i, count any emission above pi, restart the shower 
(off from pi-i) N times, and average. This will give the sum of all contributions to 

one i in eq. ( p3.39| ) 

For the O {al (hr)) —terms of the expansion of the last no-emission probability (i.e. 



the last two lines of eq. (B.39)), the function is included by only increasing the number 
of counted emissions for a trial emissions is above the merging scale. 

In our implementation, we will generate all trial emissions with running scales, and 
count relevant trial emissions with the weight 



_ Qs (pr) 
9e — / \ 9pdi,e 
as (Pe) 

Xi-lfi-l{Xi-l,pe) 



5pdf,e 



(B.43) 

if the trial emission was produced in ISR, 
if the trial emission was produced in FSR, 



with final state recoiler, 



(B.44) 



/i-l(Xi_i,pe) 



if the trial emission was produced in FSR, 
with initial state recoiler, 



where pe is the evolution scale of the trial emission, ye the energy splitting, and e the flavour 
of the initial line after the trial emission. This weight will give trial emissions generated 



^^Even if we had an analytic way of generating the integrals, implementing the average number of emis- 
sions is superior, since using trial emissions will capture correlations between potentially radiating dipoles, 
and automatically include phase space constraints. 
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with fixed scales, as can e.g. be verified for initial state splittings by using 



dp / dy as {p) 




V 



V — > — ^ — _|_ , I V X 9; where V — ' ' 



2vr xtift-ii^tvPF) V J p. P Juj y 



Weighting every trial emission with gi will exactly cancel the Qi factor in the sum, thus 
producing the desired fixed-scale terms. 

In conclusion, the algorithm to generate all O (a^ {p-r}) —terms in one multiplicity i is 

1. Start a trial shower off state S+i-i at scale pi-i. 

2. If the trial shower yields an emission, and 

1. i — 1 < n, count the emissions with weight ge if Pe > Pi 

2. i — 1 = n, count the emissions with weight ge if t {Se, Pe) > Pms- 

3. If the trial emission has been counted, restart the trial shower off state 5+i_i, with 
a starting scale pe- Repeat steps 2 and 3. 

If Pe < Pi, or Pe < Pms, or no trial emission has been constructed, set the number of 
emissions to the sum of weights, and exit. 

The average is generated by restarting this algorithm N times, and dividing the sum of 
the results by N. To generate the sum of all O [al (pr)) —terms in eq. ( |B.39| ), we step 
through the reconstructed PS history, and subtract, for each reconstructed state S'+j-i 
{i — 1 < n), the average number of emissions between and pi, as generated by the 
above algorithm. 



B.4 Summary of vi^eight generation 



This section is intended to collect the results of appendices B.l, B.2 and B.3, and to 
summarise how the necessary weights are generated. In NL'^ and UNLOPS, tree-level 
samples (Bm) are, respectively, defined as 



T' = Br 



Wm 














B, 



1 






w 














The weights Wn and w'^^ differ in that does not contain the "last" no-emission probability 
(n5+„) present in Wn- Once all approximate O (a^ (/i/j)) and O (a^ {pr)) terms in the 
weight of tree-level samples are removed, we can add NLO events, and still retain NLO 
accuracy. Collecting all O (ag (/^r)) and O {pr)) terms of appendices BJ, B^ and 
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|B.3| , we find 
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Wm 










+ 



2tt 



1 + asinRjh + ^ Qs(/iij) In f £^ 



27r ^ 



4vr V^i/^i 




pi-i 



dp / a;=^-i/^i(^,/iF) 



Pi-l 




Pt-i{y,l^F) 



dp / dyXnfni y ^Pf) r^f./a \ 

— I ±f±t ± ^Q{t{S+n+l,p) - P^ 

P J y Xnfn[Xn,PF) 

Pn+l Uj 



2tt 



(B.46) 



Pn+l 



The first terms in the expansion of the UNLOPS weight ([w^]q + [tti^]]^) do not contain 
terms proportional to 0, but are otherwise identical. The complete expression can be 
generated by using the reconstructed parton shower history of the matrix element input 
event S+n. Each step in the parton shower history corresponds to a fully reconstructed 
state S+i, and an associated production scale pi. Thus, we have enough information at step 
i to generate all O (og (pr)) —terms depending on the index i in BAt. Without specifying 
details, the method to generate the right-hand side of B.46 is 

1. Construct all PS histories for S'+„, select one history. Set v = 1 + as{pR)ki, and 
start stepping through the history at the lowest-multiplicity state S+q. Steps will be 
counted as i — 1, starting from i = 1. 

2. For each step i — 1, 

1. If i — 1 < n, increase v by the term due to the expansion of Og— I'atios, calculated 
with scale pi. 



-57- 



2. Increase v by the term due to the expansion of PDF-ratios containing the index 
i — 1. This means adding, at each step i — 1, two numerically integrated terms. 

3. Decrease v by the term due to the expansion of no-emission probabilities con- 
taining the index i — 1. This means subtracting, at each step i — 1, the average 
number of emissions between the scales pi-i and pi. In UNLOPS, only subtract 
until step i — 1 = n — 1. 

4. To arrive at (B„) samples, subtract v from the weight Wn {w'n)- Use the weight 

Wn — V [ or w'n — v) 

to fill histograms. 

In UNLOPS, there is an additional complication in that the contributions 



►n— 1 



►n— 1 



-n,n+l 



'n-s>n-l 



Br 



>n— 1 



{< - 1} 



have to be generated. For this, we generate the weights w'^ and the necessary subtractions, 
and afterwards integrate over the phase space of the n'th jet, as outlined in the generation of 



class C in section 3.2.1. Having the weights Wn, w'^, and the first terms in their expansions 

{[Wn 



[w'n]^ , at our disposal, we can produce NL or UNLOPS predictions 



for merging multiple NLO calculations by following the steps in 3.1.1 and 3.2.1 



C. Derivation of NL^ 

The aim of this section is to give a motivated derivation of the NL^ method. Since this 
derivation will explicitly require NLO-correctness as condition, the following can also be 
seen as a proof of the validity of the scheme. We will use the notation of sections H and 



A.3 



To include the parton shower resummation in a CKKW-L style, corrective weights will 
have to be applied to events. When deriving a NLO merging method, the aim must be 
that the scheme 

(a) Is correct to next-to-leading order for all exclusive n— jet observables; 

(b) Keeps the parton shower (i.e. CKKW-L) approximation for all higher orders. 

(c) Shows small dependence on the separation between matrix element and parton shower 
region, especially in the inclusive cross section. 

To handle both inclusive and exclusive NLO cross sections at the same time, let us introduce 
the symbols 

if inclusive NLO cross sections are used, 
if exclusive NLO cross sections are used. 



Bn+i^n if inclusive NLO cross sections are used, 
if exclusive NLO cross sections are used, 



(C.l) 
(C.2) 
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where B„ are the tree-level cross sections. If all the corresponding event samples are 
multiplied with corrective weights, conditions (a) and (b) read 



BnWB + L„u;l - SnWs 

Bn + Vn + In+l|n + 



PMS 

d^'rad (B 



n+l\n 



D 



oo 

n+l|n) + Bn ^ ^Wn 
1=2 



and 



Bn+lW^n = B. 



oo 

n+1 ^ ^Wn 



(C.3) 
(C.4) 

(C.5) 



i=Q 



The last equation is trivially fulfilled if we choose to reweight higher-multiplicity tree-level 
matrix elements as in CKKW-L. If we do not have control over the NLO calculation, we 
cannot assume that tree-level, virtual and real emission samples are evaluated at identical 
n-jet phase space points. Thus, without having the actual functional form of the matrix 
element weights, the merging conditions for Bn, Ln and Sn have to decouple, since we 
cannot allow wbiWl and ws to be functions of the matrix elements^^. To accommodate 
the merging constraint eq. ( |C.3| ), we make an ansatz 



i=l 
oo 

aL,o + y^^L 

i=l 

oo 

ws = asfi + ^ hs, 



WB 



WL 



iOL. 



1=1 

oo 



(C.6) 



i=l 



i=l 



We choose this form to allow for complete generality. Negative powers of ag are included 
to allow for a simple visualisation of weights stemming e.g. from division by an all-order 
expression, if such factors should be desirable. 

If we insert eq. ( |C.6 ) into eq. ( p.3| ), remember that Bn is of 0{d^), that Ln contains 
a Born term of 0{a^) and corrections of 0{a^'^'^), and that Sn is of 0(0"+-^), we can read 
off constraints on the coefficients order by order. This leads to weights of the form^^ 



WB = Wn - \wj^ ^ - 



4=1 



(C.7) 
(C.8) 
(C.9) 



^®This is not the case for MEPS@NLO in SHERPA, where full control of the matrix element functions is 
available, thus opening other avenues for NLO merging [26,27]. 

^^If L„ does not contain an additional Born term, we would not have to subtract the term [wn]^ in wb, 
which would have the benefit of fewer negative weights. 
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So far, we have allowed the coefficients c to be non-vanishing. If we naively do so, we 
allow changes of ©(a""*) to the n— jet cross section. Since the exact result (in ©(a" ) and 
0(0"+^)) should not be changed by numerically large terms, we are bound to enforce the 
conditions 



cb,i = 







ii>2) 



(C.IO) 
(C.ll) 



One way to include this condition is by replacing CB,i in the tree-level weight, which gives 
allowed weights of the form 



WB = Wn- 


Wn 
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(C.12) 



Finally, if we choose to use inclusive NLO cross sections for L„, there are non-trivial cancel- 
lations between L„ and S^, since S„ was introduced as an explicit phase space subtraction. 
If we choose wl and ws differently, these cancellations are jeopardised in higher orders. 
We thus think it reasonable to only allow the weights 



WB = Wn - [wn] ^ ~ ["^^ 1 ~ ^ \a') 

00 / 1 \* 
w, = ws = l + Y.c,, i^-J 



(C.13) 
(C.14) 



This still allows for some arbitrariness, since c^^j is not fixed. We choose a pragmatic 
approach, and exclude weights that are not easily generated by the PYTHIA8 shower^*. An 
example for such are weights with negative as order. Thus, we set 



Wb = Wn 
WL = Ws = I 



Wr 



Wr 



(C.15) 
(C.16) 



We will not reweight any ©(a""*"^)— terms. This immediately implies that the merging scale 
should be defined in the parton shower evolution variable, since otherwise, Sudakov factors 
would have to be multiplied in regions of tMs-unordered splittings [25]. Sudakov factors 
can be represented by a power series in positive powers of Qs, so that even in eq. (|C.14| ), 
we could not easily accommodate such factors. The main constraint on wl and ws is 
condition (b), which can also be interpreted as the statement that only tree- level samples 
are allowed "seeds" for higher order contributions. 

So far, we have only been concerned with conditions (a) and (b). Condition (c) be- 
comes important when combining different jet multiplicities. In NL^, this combination 
is constructed by simply summing all reweighted n— jet samples. Let M be number of 



^*We would like to point out that the approach of Platzer [32] does indeed use a smart choice to generate 
different weights, which contain Sudakov-form-factor denominators. 
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additional jets in the highest-multiplicity NLO calculation, and let us use 















= L„ = B„ 




= Sn — ~ 



B 



n.+l— ^-n. 



(n < M) 
{n > M) 
(n < M) 

(n < M) 



(C.17) 

(C.18) 
(C.19) 



The NL^ method then sums reweighted event samples: 



• If n < M, include the samples B„, L„ and S^, reweighted according to eq. ( C.15| ) or 
eq. ( p.l6|) . Remember that Sn has a negative sign. This will produce T^, and Sn- 

• If n > M, reweight B„ as in CKKW-L. This will produce T„. 
For an observable O, this produces the prediction 

n=M+l J J 



m=0 



Bm ^Wn 
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+ B„ + + I 



m ~r J-jn+l|m 
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(C.20) 
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Let us briefly investigate how this method changes the inclusive cross section in the special 
case of merging zero- and one-jet NLO calculations. At O (a^ (a^_r)), the cross section is 
given by the full NLO result by construction, while at O (og (/^r)) , we find 



:,0{S+oj) Be 
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It is useful to compare this to the result of CKKW-L 
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(C.22) 



For very low merging scales, logarithmic contributions of a single jet in B2 approaching the 
soft/collinear limit should be cancelled to better accuracy in NL^, since the one-jet NLO 
results should contain the complete logarithmic structure. However, since the zero-jet de- 
scription is given by the CKKW-L result, enhancements for the one-jet NLO contributions 
approaching the phase space boundary are not fully cancelled. For W— boson production, 
this means that NL^ does not fully compensate contributions of O In^ {pms})' 
hancements due to both jets in B2 stretching into the infrared are unchecked in CKKW-L, 
but should cancel some of the one-jet NLO terms in NL^. 

We find it difficult to assess if CKKW-L or NL'^ is more problematic. The merging 
scale value is chosen to separate the parton shower phase space from the hard matrix 
element region. Seeing that multiparton interactions at the LHC certainly play a role 
already at scales of 0(10 GeV), it is common practise to set the merging scale to a slightly 
higher value. In the particular case of W— boson production, with a merging scale of 
Pms ^ 10 GeV, double logarithms are considerably smaller than the next higher order in 
as, so that it is difficult to isolate the questionable terms. We think these important is- 
sues nevertheless, and address them, in the context of the UNLOPS method, in section |3^. 



For completeness, we add the NL^ result for exclusive NLO input, which simply does not 
contain phase space subtraction samples: 
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This ends the derivation and discussion of the NL^ method. The main conclusion of this 
section is that the allowed weights for samples in NLO merging are restricted by merging 
conditions. The constraints apply to other CKKW-L inspired NLO merging schemes as 
well. If, for example, NLO accuracy has to be safeguarded, it is mandatory to remove the 
O (ftg and O (a^ (/i/j)) -parts of the weight of tree-level events. 



-62- 



D. Derivation of UNLOPS 



In this part, we aim to give a step-by-step derivation of the UNLOPS method. We will use 
the notation defined in section and start with the UMEPS prediction for incorporating 
up to three additional jets: O as 
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Our method will be to identify which prediction for an exclusive n-jet observable we want, 
replace the UMEPS approximation by these terms, find the difference between the improved 
result and the UMEPS prediction, and remove this difference from the next-lower jet 
multiplicity. 

As a warm-up exercise, let us include a NLO calculation for zero-jet observables. We 
then want zero-jet observables to be described by the sum of Born, virtual and unresolved 
real terms, and also include the PS resummation. Keeping in mind that we do not want to 
introduce approximate 0{a1)- or 0(Qg)-terms^^, the zero-jet exclusive cross section should 
read 

Bo + 



(D.2) 
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We would thus need to remove the terms 

Bo + / Bo 



(D.3) 

from the next-lower multiplicity. Since there is no next-lower multiplicity, we simply get 
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(D.4) 



^®We use the intuitive result in eq. (CIE): If we want to include the n-jet NLO result, we have to subtract 
the 0{a'^)- and 0(Qs)-terms of the weight for n-jet tree-level events to ensure NLO accuracy. 
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'1^0 



or, in terms of the inclusive NLO cross section 
(O) = I d</.o|o(5+o,) (bo - 

+ jj 0{S+2j) (§2 - ^§3^2 

+ // / 0(5+3,) Ba I 
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Thus, we can promote UMEPS to a tree-level multi-jet merged, lowest-multiplicity NLO 
corrected calculation by simply replacing the zero-jet Born cross section with the inclusive 
zero-jet NLO cross section. Such a scheme is often called MENLOPS [38-40]. Note that 
the the total cross section is conserved, because we started from the UMEPS prediction, 
rescaled with a factor K = / Bq / / Bq. 

The derivation of a multi-jet merging scheme with NLO accuracy for any n— jet ob- 



servable is organised as follows. First, we will extend the MENLOPS result eq. (D.5) to 
simultaneously include a one-jet NLO calculation. Then, we add the two jet NLO result. 
After a short discussion, we present the general case. 

To including one-jet NLO predictions, we need to replace the UMEPS one-jet result 

by 











Bi + 


Bi 




/ B2^1 






-1,2 


. J s 



B 



3^1 
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Further we have to subtract the difference of the new one-jet prediction and the UMEPS 
case from the zero-jet part. The difference is given by 



(D.7) 
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Note that the integrated two-jet contribution 
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is integrated twice, even though the result of the first integration (iS+i) contains only 
resolved jets. Putting the pieces together, we arrive at 
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Let us rewrite this in terms of inclusive NLO calculations: 
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This result will be used in section 3.2, and is discussed there. 



Before presenting a general master formula, we will take yet another intermediate step. 



and generalise eq. (D.ll) to further include two-jet NLO predictions. Two-jet observables 
should be described by 
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To conserve unitarity, we then have to subtract the difference of the new result and the 
UMEPS contributions, i.e. 
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from the one-jet case. This gives the new one-jet contributions 
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B3^i . (D.14) 
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does not result in a state with one jet 



In the case that B2-i.i, Bg^^ or B2-> 
above the merging scale, we choose to integrate twice, as in the case of UMEPS. Thus, the 
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UNLOPS prediction for simultaneously merging zero, one and two jets at next-to-leading 
order is given by 
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Again expressing this in terms of NLO inclusive cross sections, we find 
(O) = 




(D.15) 



B 



1^0 - / B2^o 

s 



B 



1^0 



-1,2 



B 



2-s>0 



-2,3 



B 



3^0 



+ / 0{S+i,) Bi - / B2^i + 



B. 



+ / / 0{S+2j) [B2 + 

OiS+,,) §3 1 













+ 


Bi 




/ B2^l 








-1,2 


J s 


-2,3 -Is J 








) 




23 


" J 









+ 



(D.16) 















" / bLi 










Bi 


3 


/ B2^l 






/ §3^1 




3 




. J s 


3 


.Js 


3 


. J s 



Let us have a closer look at eq. ( D.15| ), and in particular how the O {a^ {hr)) —term of the 
one-jet descriptions: 



(D.17) 



The first term is the parton shower approximation of unresolved emissions in the underlying 
zero-jet configurations. The second and third terms give an approximation of unresolved 
contributions in one-jet states. Compared to the UMEPS result (J^ B2-i.i [''W2]i) ; ^^^^ should 
give an improved description. The last term in eq. ( D.17| ) is unchanged compared to 
UMEPS, and should not induce logarithmic terms in /?ms- 



Coming back to the first term in eq. ( D.17 ), it is natural to ask if reweighting the 



©(Og)— part of Bi would not give a better description. Since reweighting exclusive NLO 
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events (which contain both 0{a^) and 0(0"^^) parts) will mix higher-order terms in a 
difficult way, this interesting possibility is not examined here. We hope to come back to 
this issue for comparisons between different NLO merging prescriptions^^. 

We hope it is clear that the UNLOPS method preserves NLO accuracy by construc- 
tion, and improves the higher-order description of UMEPS further. The formal accuracy 
of exclusive n— jet observables, however, is not better than next-to-lcading order combined 
with PS resummation. Only for a limited set of observables do parton showers capture 
more than leading logarithmic enhancements. We have presented UNLOPS both for ex- 
clusive and inclusive NLO cross sections. This is motivated by trying to accommodate 
NLO calculations while requiring only minor - or ideally no - changes to the actual NLO 
implementation. 

After explicitly deriving the UNLOPS scheme for describing up to two jet observables 
next-to-leading order accuracy, we now give the UNLOPS master formula, when using 
exclusive NLO samples 
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Furthermore, since most results in this publication are produced using inclusive NLO sam- 



'We in particular think about comparisons with MEPS@NLO, were B— events are reweighted. 
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pies, we also give the UNLOPS prediction for inclusive input 
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Although the number of contributions in UNLOPS becomes somewhat unwieldy, we still 
only require M + N input event files, since some files can be reused - as in NL^. The 
procedure (including processing some input events multiple times) has been implemented 
in Pythias, and will become available in the near future. 

D.l Upgrading one-jet UNLOPS to a NNLO matching scheme 

The UNLOPS scheme has the advantage that the lowest-multiplicity cross section is not 
reweighted. This makes replacements of this term with more accurate calculations relatively 
easy. Here, we would like to hint at how UNLOPS could be shaped into a NNLO matching 
scheme. The starting point is again UMEPS, but instead of multiplying every UMEPS 
contribution with a NLO ii'-factor, we rescale with a NNLO X-factor K' . Apart from this 
change, we directly move to the UNLOPS prescription including zero- and one-jet NLO 
calculations. Then, we assume that an exclusive zero-jet NNLO calculation is available, 
producing phase space points with the weight Bq. The weight Bq should be the sum of 
the Born approximation, one-loop corrections, unresolved single real corrections, two-loop 
corrections, one-loop corrections with an additional unresolved jet, and double unresolved 
double real radiation contributions. _ 

We replace Bq in eq. ( 3.11 ) by Bq, and remove all other O (og {hr)) —terms in the 
zero-jet part: 
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We see that the inclusive cross section is given by 
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Zero-jet observables are correct to O (/^_r)), as is the description of one- and two-jet 
observables. It is of course possible to improve higher orders by including additional matrix 
elements in UMEPS-fashion. The major obstacle for implementing this method is the lack 
of available ME generators generating phase space points according to Bq. 

E. NLO merging and multiparton interactions 

Observable jets at hadron colliders are not only produced in a single energetic interactions, 
but also emerge from additional scatterings of other proton constituents. Multiparton 
interactions (MPI) models are essential in describing hadron collider data which include 
these "underlying events" [50-52]. When trying to describe the underlying event at the 
LHC, one is in practise still largely forced to use phenomenological models, although efforts 
are under way to construct more solid theoretical foundations (see [53] for a recent review) . 
Due to continuous development, it is valid to say that current phenomenological models 
offer a good description of a wide range of experimental data. 

A sophisticated MPI machinery has always been a cornerstone of PYTHIA8 [54-57]. 
Multiparton interactions in PYTHIA8 are modelled by including QCD 2 — > 2 scatterings 
in addition to the hard process. It is reasonable to assume that energetic secondary scat- 
terings induce constraints on how much beam energy would be left for further initial state 
radiation. PYTHIA8 incorporates such phase space constraints by interleaving MPI with 
parton showering: An energetic secondary scattering is produced before soft radiation. 
This is achieved by generating initial state radiation, final state radiation, and MPI in 
one decreasing sequence of evolution scales. One benefit of this method is that for a high 
shower starting scale, more jct-likc MPI are produced, which increases the Tindcrlying event 
activity for increasing hardness of the core scattering - a phenomenon called pedestal effect. 

From the point of precision QCD, we have to recognise that for observables which 
are influenced by multiple interactions, the formal accuracy of any merging method will be 
governed by MPI. Only if the influence of MPI is negligible are statements about the formal 
accuracy of the result reasonable. However, suppressing hard multiparton interactions leads 
to an inferior data description. Following the philosophy of [5] and [6], we will sacrifice the 
formal accuracy label of the NLO merging method in regions where MPI are important. 
This does not mean that we undo any improvements of our method, but only that we can 
no longer claim a particular accuracy, even in the presence of improvements. 

We include MPI in NL^ and UNLOPS in the same way as was previously done for 
UMEPS in [6] and we refer to that publication and [5] for more background. First we 
amend the normal PS no-emission probabilities with no-MPI factors, H^l'^^. This means 
that all event samples with n partons (B^, V„, I„, and L„) are multiplied with the no-MPI 
probabilities 

m—l 

nnM,7(ft,Pm), (E.l) 

i=0 

which are easily incorporated in the trial showers. Note, however, that MPI emissions 
are not taken into account when calculating the 0{al)— term of no-emission probabilities 
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(which is natural since MPI's are of 0{a^)). Then, when the shower is started from the 
reweighted (and possibly reclustered) states, using /?„ as starting scale, MPI's are included. 
As before, for n < N, any parton emission above Pms are vetoed, but if a MPI is generated, 
it is always accepted, and no emissions are vetoed in the subsequent showering. 

As in the UMEPS and CKKW-L methods, this means that any event where the n < N 
hardest jets are bove the merging scale and are from the primary interaction a, these will 
be described by the corresponding tree-level ME. In addition if n < M these jets will be 
described by the corresponding NLO ME. In both cases, the higher order Og-terms will 
be resummed to the precision of the shower. Again we note that the NLO-prediction 
will be modified by the inclusion of MPI. The modification is beyond the "leading twist" 
approximation of the NLO calculation, but may nevertheless be large, especially for jets 
with low transverse momenta. 

Before concluding, we would also point out that in this article, we have used CTEQ6M 
PDFs in the generation of secondary scatterings. This is not advisable, since the domi- 
nant contribution to the underlying event stem from soft secondary scatterings. For such 
scatterings, PDFs are evaluated at low scales 0(1 GeV), and very small x— values, i.e. in 
a region where NLO PDFs are poorly constrained and need not even be positive definite, 
which clearly is problematic in the probabilistic MPI picture. When developing a future 
tune to be used together with NLO merged predictions, we will utilise NLO parton dis- 
tributions for the hard interaction, while employing leading-order PDFs in the multiple 
interactions and parton showers. 
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